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ABSTRACT 


This  report  describes  work  undertaken  to  resolve  the 
SNAP-50/SPuR  iiquid-meta!  boiler  two-phase  flow  stability  problems. 

A  review  is  presented  of  existing  literature  pertaining  to  forced- 
convection  vaporization  heat  transfer  and  two-phase  flow  stability. 
Dimensionless  simulation  criteria  are  defined  on  the  basis  of  the 
significant  parameters  determined  from  the  survey.  The  test  pro¬ 
gram,  in  whfCh'  freon  113  was"used  as  the  test  fluid  and  the  values 
of  the  dimensionless  two-phase  flow  stability  parameters  applicable 
to  the  SNAP'50/SPUR  boiler  were  simulated,  is  presented.  The  design 
of  the  transparent  closed  test  loop,  which  incorporates  a  five- tube 
boiler,  is  described.  Conciusions  and  discussions  are  presented, 
based  on  the  measured  data  and  the  high  speed  photographs  of  two- 
phase  flow  in  both  plain  tubes  and  tubes  with  twisted  tape  inserts, 
which  indicate  that  stable  flow  can  be  achieved  in  the  SNAP-50/SPLIR 
boiler  by  inserting  at  the  entrance  of  each  boi ier  tube  an  orifice 
which  has  a  pressure  drop  between  100  and  200  percent  of  the  two- 
phase  pressure  drop. 


PubPcat'on  of  this  technical  documentary  report  does  not 
constitute  Air  Fo'ce  app'-oval  of  the  report's  findings  or  conclusions. 
It  is  published  only  *or  the  exchange  and  st  i.mu  lat  ion  of  ideas. 
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NOMENCLATURE 


Cp  Specific  heat,  Btu/lb-°P 

0  Tube  diameter,  ft 

f  Frequency,  cps 
G  Mass  velocity,  Ib/hr-sq  ft 

Heat  transfer  coefficient,  Btu/hr-sq  ft-®F 
^fg  vaporization,  Btu/lb 

L  Tube  length,  ft 

F  Pressure,  psi 

Q/A  Heat  flux,  Btu/hr-sq  ft 
P  Density,  Ib/cu  ft 

T  Temperature,  ®F 

U  Overall  coefficient  of  heat  transfer,  Btu/hr-sq  ft-®F 

W  Flow  rate,  Ib/hr 

^  Qua  1 i ty 

Subset'  ipts 


boi ler 


c  critical 


^  freon 


in  in 


i  iquid 


out  out 


sat  saturation 
TP  two- phase 
tube  tube 


vapor 
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SECTION  I 


INTRODUCTION 


The  SNAP-50/SPUR  power-conversion  system  requires  a  boiler  that  converts 
subcooled  liquid  potassium  Into  dry  potassium  vapor  using  liquid  lithium, 
heated  in  a  nuclear  reactor,  as  the  heat  source.  The  boiler  is  a  shell-and- 
tube  design,  with  the  potassium  vaporizing  inside  parallel  tubes  which  contain 
twisted  tape  Inserts.  This  report  describes  work  undertaken  to  resolve  the 
SNAP-50/SPUR  liquid-metal  boiler  two-phase  flow  problems  by  simulating  flow 
conditions  with  conventional  fluids.  Two-phase  flow  problems  may  arise  due 
to  flow  stability  phenomena,  multichannel  flow  distribution  effects,  or  high 
quality  dry-wall  transition  effects,  as  des>,'ibed  below. 

•  Flow  Stabi 1 i ty.  Flow  stability  is  a  problem  inherent  in  forced- 
convection  boilers  of  the  SNAP-50/SPUR  type,  in  which  subcooled  liq¬ 
uid  is  passed  into  parallel  passages  and  vaporized  to  high  quality. 
Small  local  changes  in  pressure  drop  or  heat  transfer  rate  can  cause 
rapid  and  substantial  shifts  in  flow  rate  through  individual  tubes. 
One  reason  for  this  is  that,  in  the  two-phase  region,  the  curve  of 
pressure  drop  vs  flow  rate  can  have  a  negative  slope  owing  to  the 
phase  change.  That  is,  a  decrease  in  flow  leads  to  an  increase  in 
vapor  quality  and,  since  the  pressure  drop  associated  v/ith  vapor 
flow  is  much  greater  than  that  for  liquid,  a  greater  pressure  drop 
is  required.  The  imposed  pressure  drop,  however,  tends  to  remain 
constant;  therefore,  a  decrease  in  flow,  which  requires  an  increased 
pressure  drop,  in  turn  leads  to  a  further  decrease  in  flow. 

•  Multichannel  Flow  Distribution.  Flow  distribution  is  another  prob¬ 
lem  existing  in  any  heat  exchanger  with  multichannels.  Poor  flow 
distribution  causes  different  rates  of  temperature  or  quality  in¬ 
crease  in  each  tube,  so  that  the  required  exit  condition  is  not 
attained  at  an  equal  distance  from  the  inlet  in  every  tube.  The 
effect  of  unequal  distribution  is  to  require  that  every  tube  be 
lengthened  so  that  each  tube  vm  1 1  produce  a  satisfactory  exit  con¬ 
dition;  this  increases  the  boiler  weight.  Good  flow  distribution 
normally  is  promoted  by  the  design  of  the  inlet  iranifold. 

•  High  Quality  Dry-Wall  Transition.  At  high  quality  in  forced  convec¬ 
tion  vaporization,  there  is  a  transition  from  a  liquid  covered  v/al  1 
to  a  partially  dry  wall,  with  a  corresponding  decrease  in  heat  trans¬ 
fer  rate.  The  quality  at  which  this  occurs  cannot  be  predicted, 

and  the  factors  affecting  this  transition  and  the  flow  patterns  are 
uncert'-in.  It  is  not  known,  for  example,  v;hecher  the  small  volume 
of  liquid  present  covers  one  percent  of  the  wall,  or  none  of  the 
wall.  The  use  of  whirlers  to  centrifuge  the  liquid  to  the  wall  may 
increase  the  amount  of  wall  covered  by  liquid.  Observation  of  the 
flow  pattern  for  various  whirlers  with  high  speed  motion  pictures 
is  a  very  valuable  tool  for  clarifying  this  important  area.  Essen¬ 
tially  every  advancement  in  the  field  of  boiling  heat  transfer  and 
two-phase  flow  has  been  the  result  of  visual  observations. 


Manuscript  released  by  author  October  1966  for  publication  as  an  RTO  Technical 
Documentary  Report. 


Since  liquid  the''ma!  conduct  •  v  i  t /,  which  is  the  only  p'-ope’'ty  of  liquid 
metals  that  varies  greatly  f''om  the  va^ue  fof  common  fluids,  does  not  directly 
influence  the  two-phase  flow  pi'ocesses,  't  is  possible  to  simulate  liquid- 
metal  flow  processes  with  corventiona*  flu'ds.  “his  approach  allows  the  use 
of  simple  flow  visual ization  techniques,  which  provide  a  mo^e  complete  under¬ 
standing  of  the  majo'  p'-oblem  areas,  and  cons  i  de'ab!  y  reduces  test  complexity 
during  initial  evaluation  and  solution  of  the  two-phase  flow  problem  areas. 

The  objective  of  this  Investigation  v/as  to  proceed  as  far  as  possible 
toward  the  solution  of  the  SNAP  50/SPUR  boiler  two-phase  flow  problems  using 
information  available  in  the  literature  and  conventio..a  ■  fluid  testing.  This 
approach  minimized  the  amount  of  liquid-metal  testing  required  to  obtain  the 
final  solution  to  these  problems,  and  thus  'esulted  in  a  significant  cost 
sav  i  nc. 

This  report  contains  a  b'ief  Ifte'atu'e  'eview  that  produced  information 
on  the  important  parameters  which  we'e  then  used  to  formulate  flow  simulation 
.riteria.  The  test  program  resjjting  T^om  the  application  of  these  simulation 
criteria  and  the  transparent  test  loop  designed  to  measure  the  effect  of  che 
controlling  parameters  are  described  finally,  the  results  of  the  test 
program  are  presented  and  discussed. 
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SECTION  2 


TECHNICAL  DISCUSSION 


GENERAL 


Flow  instability  may  occur  under  an  unkr.own  combination  of  heat  exchanger 
boundary  conditions,  flow  rate  vs  pressure  drop  relationships  (functions  of 
fluid  properties),  and  external  disturbances,  and  has  been  detected  by  many 
investigators.  It  is  clear  from  the  experimental  evidence  that  destabli;?ed 
flow  conditions,  both  with  and  without  heat  addition.  For  example,  Prandt I 
showed  that,  in  the  transition  region,  single-phase  flow  in  a  pipe  could 
destabilize  and  oscillate  between  laminar  and  turbulent  flow.  Of  greater 
concern  In  this  program,  however,  are  the  flow  instabilities  associated  with 
two-phase  flow;  only  these  cases  will  be  discussed  in  succeeding  sections. 

Understanding  and  analysis  of  the  flow  instability  problem  requires  a 
knowledge  of  many  related  technical  specialties.  Some  of  the  instabilities 
are  the  result  of  flow  or  heat  transfer  regime  transitions,  and  nvJny  can  be 
qualitatively  explained  by  more  than  one  mechanism.  A  thorough  understanding 
of  all  heat  transfer  mechanisms  and  flow  regimes  possible  with  forced-convection 
vaporization  heat  transfer  and  two-phase  flow  is  required  to  preclude  erroneous 
conclusions.  Therefore,  a  brief  discussion  of  the  state  of  the  art  in  forced- 
convection  vaporization  heat  transfer.  Including  research  results  from  an 
earlier  phase  of  this  investigation,  is  presented  below,  prior  to  a  review  of 
the  literature  containing  Information  on  two-phase  flow  stability. 

FORCED-CONVECTION  VAPORIZATION  HEAT  TRANSFER 


A  large  number  of  both  experimental  and  analytical  investigations  of  forced- 
convection  vaporization  heat  transfer  have  been  performed  during  the  past  20  years. 
Reference  I*  contains  an  excellent  review  of  the  state  of  the  art  up  to  1957.  The 
process  depends  upon  a  great  number  of  variables  and  is  made  more  difficult  by 
the  canplexity  of  the  various  two-phase  flow  patterns  which  occur  as  the  quality 
increases  along  the  tube. 

Because  of  the  large  differences  in  the  various  flow  patterns  which  occur 
as  the  volume  fraction  of  the  vapor  increases  along  the  tube,  a  separate  cor¬ 
relation  is  required  for  each  flow  pattern.  Therefore,  familiarity  with  the 
characteristics  of  each  flow  pattern  is  required  before  an  analytical  model 
can  be  formulated  to  yield  correlations  for  the  heat  transfer  coefficient  In 
each  regime,  and  the  quality  at  which  flow  regime  transitions  occur. 


Visual  sjdies  of  adiabatic  two-phase,  two-component  flow  have 
provided  information  on  the  flow  patterns.  Attempts  have  been  made 
develop  flow  regime  maps,  which  would  make  it  possible  to  predict 
the  boundaries  between  the  observed  flow  patterns  While  the 


to 


♦Numbers  refer  to  references  at  the  end  of  this  report. 
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flow  regimes  depend  to  some  extent  upon  whether  the  flow  is  vertical  or 
horizontal,  the  flow  regimes  shown  in  Figure  I  appear  to  Lc  nost  common  for 
both  orientations.  The  flow  patterns  shown  are  as  follows; 

a.  Bubbly  Flow- -This  pattern  occurs  at  very  low  quality  and  consists  of 
individual  bubbles  of  vapor  entrained  in  the  liquid  flow.  The 
bubbles  are  formed  at  nudeation  sites  on  the  wall. 

b.  Plug  or  Slug  Flow--As  the  qualify  increases,  the  individual  bubbles 
agglomerate  and  grow  to  form  plugs  or  slugs  of  vapor,  which 
periodically  pass  a  given  point  on  the  wall. 

c.  Annular  Flow — At  higher  quality,  the  vapor  flows  in  a  high-velocity 
core  in  the  cente-  of  the  pipe,  and  a  thin  film  of  liquid  covers  the 
pipe  wall.  Liquid  droplets  are  generally  dispersed  in  the  vapor 
core. 

d.  Mist  Flow--At  moderate  to  high  quality,  the  liquid  film  on  the  wall 
disappears  and  all  the  liquid  is  dispersed  throughout  the  vapor  as 
droplets. 

Forced-convection  vaporization  heat  transfer  data  have  shown  that  the 
heat  transfer  coefficient  for  a  fixed  geometry  and  flow  rate  increases  with 
increasing  quality^^>^X  At  a  critical  quality,  which  varied  considerably  with 
fluid,  geometry,  and  heat  flux,  there  is  a  sharp  decrease  In  the  heat  transfer 
coefficient  toward  the  value  characteristic  of  vapor  flow,  as  shown  in 
Figure  2.  When  the  quality  reaches  100  percent,  the  coefficient  equals  the 
vapor  heat  transfer  coefficient.  Numerous  heat  transfer  correlations  have 
been  recommended  vl>^)6,7)  for  various  flow  regimes. 

It  is  often  of  great  importance  to  be  able  to  predict  the  quality  at 
which  the  sharp  decrease  in  heat  transfer  coefficient  occurs,  for,  if  the  heat 
flux  is  the  imposed  condition,  burnout  may  occur  beyond  this  point.  It  has 
been  suggested'^/  that  this  quality  corresponds  to  the  point  at  which  the 
transition  from  annular  flow  to  mist  flow  occurs.  This  is  very  likely,  since 
this  is  a  transition  from  a  liquid  covered  wall  to  a  gas  blanketed  wall,  which 
is  similar  to  the  pool  boiling  burnout  condition  that  occurs  in  changing  from 
nucleate  to  film  boiling.  Many  reasons  have  been  suggested  for  the  transition 
f'om  annular  to  mist  flow,  and  numerous  correlations  of  the  forced-convection 
vaporization  burnout  condition  have  been  attempted. 

Most  of  the  visual  studies  of  forced-convection  vaporization  have  been 
performed  under  adiabatic  conditions  because  of  the  relative  ease  of  this  type 
of  experiment.  In  addition,  most  of  the  heat  transfer  tests  have  used 
electrical  heating,  thereby  imposing  the  heat  flux  boundary  condition.  An 
experiment^®)  performed  at  AIResearch  during  an  earlier  phase  of  this  program 
used  high  speed  motion  pictures  to  study  the  forced- con vect ion  vaporization 
process,  when  the  overall  temperature  difference  was  the  imposed  boundary 
condition.  Motion  pictures  were  obtained  of  all  the  significant  flow  patterns; 
the  resulting  film  may  be  borrowed  from  the  writer. 
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Figure  2.  Cliai  acter  i  f.  t  ic  F'^rced-Convcct  ion  Vaporization 
Heat  Transfer  Coefficient 
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The  following  conclusions  were  bosed  on  resuits  of  an  earlier  phase 
of  this  program,  which  Is  reported  in  Reference  8  (included  as  an  appendix 
to  this  report). 

Bubbly  Flow 

Because  of  the  high  ratio  of  liquid  to  vapor  density  for  most  fluids  of 
interest,  the  vapor  mass  fraction  (quality)  in  the  bubbly  flow  regime  is 
generally  much  less  than  I  percent.  The  bubbly  flow  regime  can  occur  in  a 
subcooled  liquid,  since  the  bubbles  nucleating  at  the  wall  may  depart  from  the 
v/a  1 1  because  of  the  inertia  caused  by  their  growth  rates.  Thermal  equilibrium 
does  not  exist,  as  there  is  insufficient  time  for  the  bubbles  to  condense 
before  the  surrounding  liquid  is  heated  to  the  saturation  temperature  as  it 
flows  along  th®  tube. 

In  the  bubbly  flow  regime,  the  volume  fraction  of  the  vapor  is  very  small, 
so  that  it  has  little  effect  on  the  liquid  velocity.  Bubbles  grow  and  depart 
from  nucleation  sites  on  the  wall  much  as  they  do  in  nucleate  pool  boiling. 
Depending  upon  the  wall  and  fluid  temperatures  and  the  nucleation  character¬ 
istics  of  the  wall,  the  nucleate  pool  boiling  heat  transfer  mechanism  may  or 
may  not  dominate  the  forced-convection  heat  transfer  to  the  liquid  occurring 
between  the  bubble  nucleation  sites  In  the  very  narrow  quality  range  over 
which  bubbly  flow  exists,  it  seems  reasonable  to  assume  that  the  heat  transfer 
can  be  predicted  by  superimposing  liquid  forced  convection  and  nucleate  pool 
boiling,  provided  the  superheat  temperature  difference  is  not  great  enough  to 
produce  film  boiling. 

Plug  Flow 

As  the  vapor  volume  fraction  or  the  tube  length-to-d iameter  ratio 
increases,  the  individual  bubbles  grow  and  agglomerate  to  form  plugs  or  slugs 
of  vapor,  which  occur  periodically  along  the  tube.  Even  though  the  vapor 
volume  fraction  has  increased,  the  corresponding  quality  is  still  very  low. 
However,  the  presence  of  the  larger  volume  fraction  of  vapor  begins  to  cause 
a  significant  increase  in  the  liquid  velocity  in  the  plug  flow  regime. 

Bubbles  may  nucleate  and  depart  from  the  wall  at  numerous  sites  inside 
the  tube  in  the  plug  flow  regime.  The  heat  transfer  mechanism  is  probably 
the  same  as  that  in  the  bubbly  flow  regime,  a  superposition  of  forced 
convection  to  a  liquid  and  nucleate  pool  boiling. 

The  plugs  of  vapor  are  compressible  volumes  which  make  possible  flow 
oscillations  within  the  tube  even  though  the  flow  entering  the  tube  is  steady. 
Actual  flow  reverfial  was  observed  in  the  motion  picture  study.  It  is  possible 
that  the  periodic  nature  of  plug  flow,  combined  with  the  campressible  volume 
it  provides,  is  the  source  of  many  of  the  exit  flow  oscillations  that  occur 
in  forced-convection  vaporization, inasmuch  as  all  other  flow  regimes  are  steady. 

Annular  Flow 

Theoretical  stability  analysis  of  a  liquid  film^*^^  predicts  tnat  a  liqui 
film  is  stable  for  sufficiently  small  Reynolds  numbers,  independent  of  the 
vapor  velocity.  Since  the  liquid-film  Reynolds  number  decreases  as  the 
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Mua'ity  increases  in  a  forced-convection  vaporizer,  the  liquid  annulus  will  be 
stable  at  sufficiently  high  quality,  independent  of  the  value  of  the  vapor 
veloc i ty. 

While  annular  two-phase  flow  is  sometimes  idealized  so  that  all  the  liquid 
is  in  the  annular  film  and  all  the  vapor  in  the  high  velocity  core,  observa¬ 
tions  generally  have  shown  that  droplets  of  liquid  exist  in  the  vapor  core. 

This  is  consistent  with  the  stability  analysis  which  predicts  that  a  high- 
Reynolds-riumber  liquid  film,  such  as  exists  at  low  quality,  may  be  unstable 
and  produce  droplets.  These  droplets  may  or  may  not  rejoin  the  liquid  annulus, 
depending  upon  the  magnitude  of  the  turbulent  vapor  fluctuations,  the  transit 
time  of  the  droplets  through  the  vaporizer  and  the  passage  length-to-diameter 
ratio. 

The  liquid  vapor  interface  between  the  continuous  liquid  film  on  the  wall 
and  the  high-velocity  vapor  core  is  wavy.  Although  there  are  active  bubble 
nucleation  sites,  these  probably  have  .j  negligible  effect  on  the  rate  of  vapor 
generation.  The  vapor  is  generated  primarily  by  vaporization  from  the  liquid 
vapor-core  interface  inside  the  tube,  and  not  by  the  formation  of  bubbles  at 
the  wall.  The  heat  transfer  in  annular  flow  is  very  similar  to  that  for 
condensation  inside  a  tube;  heat  is  transferred  through  the  liquid  film  to 
the  interface  with  the  vapor  core. 

Annular-Mist  Flow  Regime  Transition 

A  description  of  the  transition  from  annular  to  mist  flow  is  of  great 
interest,  since  this  is  presumably  the  point  at  which  the  vaporizer  heat 
transfer  coefficient  experiences  a  sharp  decrease  and  causes  the  forced- 
convection  vaporization  burnout  condition. 

The  following  significant  observations  of  the  wall  drying  process  were 
made.  A  small  dry  spot  forms  and  grows  in  all  directions,  as  the  liquid 
vaporizes,  because  of  conduction  heat  transfer  through  the  liquid.  The  small 
strips  of  liquid  on  the  wall  remain  almost  stationary  relative  to  the  high- 
velocity  vapor  and  the  liquid  droplets  flovying  in  the  tube.  The  dominant  heat 
transfer  mechanism  is  conduction  through  the  liquid.  Nucleation  may  produce 
the  initial  dry  spot  on  the  wall,  but  has  a  small  effect  on  heat  transfer. 

This  is  the  same  drying  process  that  would  occur  with  a  thin  film  of  liquid 
in  a  hot  pan  whose  temperature  is  not  great  enough  to  cause  nucleate  boiling. 

As  a  result  of  the  visual  observations  of  the  transition  process,  it  was 
concluded  that  the  transition  is  due  simply  to  vaporizing  of  the  liquid  film 
on  the  wall  as  a  result  of  conduction  heat  transfer  through  the  liquid,  while 
the  vapor  core  still  contains  a  dispersion  of  liquid  droplets  formed  upstream 
at  lower  quality  and  carried  along  by  the  vapor.  While  the  vol ume-f raction  of 
these  droplets  is  very  small,  they  account  for  a  substantial  mass  fraction 
because  of  the  high  1 iquid-to-vapor  density  ratio. 

In  summary,  the  liquid  droplets  dispersed  throughout  the  vapor  core  are 
produced  as  a  result  of  a  film  instability  at  low  quality.  The  transition  from 
annular  to  mist  flow  occurs  when  the  liquid  annulus  vaporizes  as  a  result  of 
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conduction  through  the  thin  liquid  film,  while  1  .,uid  droplets  remain 
dispersed  throughout  the  vapor  core.  The  quality  at  the  transition  is  deter¬ 
mined  by  the  mass  of  liquid  contained  in  the  droplets  when  the  transition  occurs. 

Mist  Flow 


In  mist  flow,  the  wall  is  dry  except  when  the  liquid  droplets  in  the 
por  core  impinge  upon  the  wall  and  vaporize.  Most  cf  the  heat  transfer  is 
from  the  wall  to  the  vapor.  After  -he  heat  is  transferred  to  the  vapor,  it  is 
then  transferred  to  the  liquid  droplets;  the  mist  flow  vaporization  process 
actually  takes  place  in  the  interior  of  the  pipe  and  not  at  the  wall.  For 
this  reason,  the  temperature  of  the  vapor  in  the  mict  flow  regin.e  is  g'eater 
than  the  saturation  temperature;  thermal  equilibrium  does  not  exist  in  the 
tube.  The  mechanism  of  heat  transfer  in  mist  fiov/  is.  therefore,  forced 
convection  to  a  gas. 

Flow  Stab i 1 i ty 

For  low  values  of  pressure  drop  across  the  throttle  valve  at  the  entrance 
to  the  heated  tube,  flov;  instability  was  observed.  The  flov/  rate  fluctuations 
apparently  were  due  to  the  interaction  of  the  two-phase  flow  pressure  drop 
characteristics  with  the  compressible  volumes  in  the  system^*^/.  At  all  heat 
fluxes  and  flov/  rates  investigated,  the  inlet  flov/  rate  could  be  made  steady 
with  a  sufficiently  great  pressure  drop  across  the  inlet  throttle  valve.  Even 
though  there  was  steady  flow  into  the  tube,  the  flow  rate  and  temperature  at 
the  exit  were  seen  to  fluctuate.  This  apparently  was  due  to  the  periodic 
nature  of  plug  flow,  the  compressible  volume  of  the  plugs,  and  the  fact  that 
thermal  equilibrium  does  not  exist  in  mist  flow. 

TWO-PHASE  FLOW  STABILITY  LITERATURE  REVIEW 

Causes  of  Flow  Instability 

The  relative  influence  of  a  given  design  parameter  on  the  stabi  I  i ty  of 
fluid  flov/  systems  depends  quite  strongly  on  the  particular  type  of  instability 
to  which  a  system  is  subject.  Nonetheless,  where  catastrophic  power  excursions 
are  not  involved,  the  literature  dealing  with  two-phase  fiow  instability 
indicates  that  the  follov/ing  parameters  are  i  nf  1  uen  t 'a  1  ;  (l  fluid  inlet  sub- 
coo  1 ing,  (b)  heat  f 1  ux,  (c )  pressure,  (d )  inlet  velocity,  (e  orifices,  and 
(f  )  acceleration  fields- 

! .  Subcool i nq 

A  primary  cause  of  instability  in  boiling  water  reactors  has  been  inlet 
subcooling,  and  this  has  been  particularly  evident  where  large  void  fractions 
v^ere  present.  Beckjord(**^  attributed  the  flov/  instabilities  he  observed  in 
a  boiling  water  reactor  to  inlet  subcooling.  It  has  also  been  noted  that, 
in  nonreactor  circulation  tests  v/i  th  heat  addition,  the  condition  for  the 
onset  of  oscillation  is  strongly  dependent  upon  the  inlet  temperature,  and 
that  the  observed  instabilities  increase  in  severity  with  increases  in  the 
amount  of  subcooling.  This  is  in  agreement  with  formal  calculations  of 
density  changes  in  coolants,  which  lead  to  the  conclusion  that  '^ubcooling 
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increases  unstable  density  changes,  it  aiso  is  In  agreement  with  other 
calculations  which  demonstrate  that  increasing  subcooling  enlarges  regions  of 
operation  in  which  the  slope 
negative,  and  hence  unstable 


of  the  curve  for  pressure  drop  vs  flow  rate  is 
(I?. 


On  the  other  hand,  in  tests  conducted  of  boiling  burnout  with  water  In 
vortex  or  swirl  flow  ('3)^  where  degree  of  subcooling  was  treated  as  a  variable, 
excellent  stability  was  exhibited  up  to  the  burnout  point.  This  would  suggest 
that  subcooling  perhaps  does  not  influence  the  stability  of  swirl  flow,  or  that 
the  Influence  is  extremely  weak;  more  definitive  investigations  are  required 
to  confirm  or  refute  this  conjecture. 


Recent  results  of  Gouse^*^)  from  a  forced-flow,  paral lei -channel  system 
showed  that  there  was  a  range  of  subcooling  within  which  the  flow  would 
oscillate.  If  the  inlet  subcooling  was  either  greater  or  less  than  this  range, 
the  flow  was  steady. 


2.  Heat  Flux 


There  is  considerable  disagreement  regarding  the  effect  of  heat  flux  on 
stability.  This  probably  is  because  heat  flux  should  not  be  considered  an 
independent  variablej  It  can  be  determined  from  inlet  velocity,  subcooling, 
exit  quality,  and  geometry.  Some  experimental  workers  find  that  instabilities 
have  a  tendency  to  be  strongly  damped  and  reduced  in  amplitude  when  design 
heat  inputs  are  raised  sufficiently.  Although  not  clearly  understood,  the 
phenomena  seem  to  be  connected  with  the  possibility  that  damping  forces  may 
increase  more  rapidly  than  forces  driving  the  instabilities.  On  the  other 
hand,  Gouse^*^/,  Quandtv*^)^  and  Meyer(l6)  have  observed  that  oscillatory 
instabilities  can  be  induced  by  increasing  the  heat  flux. 

Investigators  apparently  have  not  considered  the  effect  on  stability  of 
a  constant  wall  temperature  boundary  condition  instead  of  a  uniform  wall 
heat  flux. 


3.  Pressure 


?  I 


?  I 


Oscillations  are  strongly  damped  and  relatively  low  in  amplitude  at 
elevated  pressures.  Blubaugh  and  Quandt'*^'  noted  that  conditions  for  the 
onset  of  oscillation  were  primarily  a  function  of  the  channel  inlet  temperature 
and  the  pressure.  Quandt's^*^)  theoretical  analysis  showed  that  flow  oscil¬ 
lations  become  less  prevalent  at  higher  pressures.  Mendler,  et  al^*®'  also 
observed  that  flow  fluctuations  in  natural  circulation  tests  increased  with 
the  reduction  of  pressure  and  inlet  temperature. 


While  this  subject  requires  more  study,  the  influence  of  pressure  seems 
to  be  connected  with  the  fact  that,  at  higher  pressures,  the  ratio  of  liquid 
to  vapor  densities  is  reduced  and,  in  a  sense,  the  two  phases  are  more  similar. 
StenningO^)  has  shown  in  his  analysis  that  oscillatory  flow  will  exist  for  a 
variety  of  configurations,  provided  that  the  density  ratio  across  the  system 
exceeds  a  critical  value, which  depends  on  the  geometry  and  the  flow  relation¬ 
ships.  Calculations  performed  by  Wright'^®^  also  indicate  that  flow  stability 
exists  if  the  specific  volume  ratio  is  below  a  critical  value. 
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Inlet  Velocity 

Inlet  velocity,  mass  velocity,  and  flow  rate,  which  are  discussed  bv 
va'^ious  investigators,  a-e  different  names  for  the  same  variable. and  all  will 
be  nrouped  under  the  above  heading.  Approximate  calculations  indicate  that 
an  increase  in  the  design  inlet  velocitv  has  a  stabilizing  influence  in 
le,  once-through  high  qua’ it/  boilers.  Jones  (21)  has  shown  that  the  tendency 

L  towa’'d  instability  varies  inversely  with  the  square  of  the  subcooled  fluid 

velocity.  Gouae'vl^)  observed  that  increasing  the  flow  rate,  with  other 
variables  held  constant,  tends  to  stabilize  the  system.  At  the  present  time, 
however,  the  sparseness  of  the  literature  does  not  lend  itself  to  definitive 
statements. 

5 .  O'^ifice  Pressure  Ofou 

Severa'  i  nvesc  i  oato  •  s  have  considered  the  effects  of  orj^^ices  at  the 
•  nlet  on  flow  stability.  f’''aasll2j  has  shown  that,  i  ■'  lov/  pressu-e  s^bcoo’ed 
wate*",  orif  icing  has  a  p-onounced  stab'Mzing  'nfluence;  moreover,  the  effect 
•ncreases  w  th  the  'atio  of  the  orifice  pressjrg  drop  to  the  pressure  drop 
fo'  the  -'est  of  the  passage.  Stenning  19)  considered  analytically  the  time- 
varying  flow  of  a  boiling  liquid  through  a  series  of  ducts  and  heaters  for  the 
case  in  wh I c h  the  flow  'ate  is  controlled  by  a  downstream  orifice.  He  showed 
.  that  oscillatory  flow  will  exist,  provided  that  the  density  ratio  across  the 

system  exceeds  a  certain  c'itical  value.  He  concluded  that  upstream  orificing 
will  help  •■o  damp  out  these  oscillations. 

A  ore-pass  boiler  a  *  so  has  been  studied  by  Jones  21)  for  the  case  in 
which  the'e  is  significant  subcooling.  He  found  that,  while  an  orifice  can 
be  a  stabiliz'ng  Influence,  it  also  can  destabilize  the  system.  From  his 
conclusions  r(;-g3rding  the  stabilizing  influence  of  high  subcooled  velocities, 

't  s  by  no  means  clear  that  orificing  would  be  effective  in  high  performance 
systems  whe  e  there  is  tu'bulent  flow  and  vigorous  subcooled  boiling  with  a 
relatively  high  friction  factor. 

(22) 

Wall'S  and  Heasley  have  discussed  a  natural  circulation  loop  with  a 
corst'ictior  at  the  top  of  the  heated  section.  They  assert  that  their  results 
have  wide  application  because  many  systems  without  outlet  orifices  have  most 
o*^  the  pressure  drop  at  the  exit  where  velocities  and  shear  stresses  are  high, 
e  ^hey  obsc-rved  instability  in  a  pentane  loop  when  the  heater  exit  was 

const' icted. 

Several  investigators  have  demonstrated  that  the  stability  of  coolant  flow 
du'ing  boil  ng  is  influenced  by  the  pump  and  system  flow  characteristics. 
Stability  's  improved  as  the  test~section  pressure  drop  becomes  a  smaller  part 
of  the  tota •  system  pressure  drop.  The  change  in  pressure  drop  caused  by  a 
pr'ase  charge  oT  the  coolant,  then,  has  less  influence  on  the  flow  rate,  since 
this  p'essj-e  drop  is  a  smaller  part  of  the  total  pump  head. 

(23 ) 

LowdermilK  found  that  it  always  was  possible  to  stabilize  the  flow  in 
a  single  tube  with  a  great  enough  inlet-restriction  pressure  drop.  The  required 
pr^ssu'e  drop  was  concluded  to  be  a  function  of  L/D  ratio  and  mass  velocity. 
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6. 


Acceleration  Field 


Depending  upon  the  direction  of  the  field,  acceleration  can  be  extremely 
important,  either  in  stabilizing  or  in  destabilizing  a  fluid  flow  system.  At 
one  extreme  is  the  well  known  catastrophic  Taylor  instability,  resulting  when 
the  field  is  in  the  "wrong  direction"  across  the  interface  between  liquid  or 
gaseous  media  of  different  densities  (24),  At  the  other  extreme  is  the  case 
of  the  highly  stable  interface  that  exists  in  high-velocity  separated  two- 
phase  flow  in  a  tight  spiral.  Because  of  the  extreme  difficulties  involved 
in  carrying  out  both  theoretical  and  experimental  work  in  reduced  gravity  fields, 
very  little  work  has  been  reported  in  this  area  beyond  the  fact  that  the 
tendencies  toward  instability  seem  to  increase  as  the  conditions  for  free  fall 
are  approached. 

Classes  of  Flow  InstabiU  r i es 


Although  authors  have  not  formally  categorized  the  hydrodynamic  instabil¬ 
ities  in  two-phase  fluid  flow  systems,  those  not  involving  actual  burnout  may 
be  divided  approximately  into  the  following  four  classes: 

•  Flow  regime  transitions 

•  Negative  derivative  of  pressure  drop  with  respect  to  flow  rate 

•  Compressibility 

•  Parallel  channel  coupling 

A  fifth,  the  amplification  of  a  "maldistribution,"  has  been  reported  by 
Adorni,  et  but  with  neither  a  description  nor  discussion  adequate  for 

review  at  this  time.  Also,  Quandt (15),  on  the  basis  of  an  analytical  and 
experimental  investigation  of  f low  osci 1 lations,  concludes  that  another  type 
of  flow  instability,  "stable  flow  oscillations,"  associated  with  definite 
frequency  and  reproducibility  from  cycle  to  cycle,  should  be  recognized. 

Systems  subject  to  different  classes  of  instabilities  may  exhibit  quite 
similar  performance.  For  example,  oscillations  may  be  associated  v/ith  flow 
regime  transitions  and  parallel  channel  coupling  as  well  as  with  compressi¬ 
bility.  Flow  excursions  may  be  associated  with  all  classes.  In  spite  of 
these  similarities,  it  is  necessary  to  treat  each  class  separately  because  of 
basic  differences  in  the  phenomena  believed  to  be  responsible  in  each  case. 

In  addition,  the  type  of  instability  Is  a  function  of  the  imposed  boundary 
conditions  on  the  heat  exchanger.  Typical  boundary  conditions  are  constant 
pressure  drop,  constant  inlet  pressure,  and  constant  Inlet  flow  rate.  The 
boundary  condition  selected  in  the  analysis  or  imposed  in  the  experiment  will 
determine  the  different  types  of  instabilities  that  will  be  possible. 

I .  Flow  Regime  Transitions 

The  recognized  types  of  flow  regime  transitions  that  can  give  rise  to 
unstable  operation  are  nucleate-to-f i Im  boiling,  bubble  flow-to-slug  flow- 
to-annular  flow  and  laminar-to-turbulent  flow.  These  transitions  give  rise 
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to  such  instabilities  as  pressu»'e  variations,  transient  pressure  drop  (and 
'•ecove''y),  and  sustained  oscillations. 

/  2a  \ 

Isbin,  et  aM  has  noted  that  large  pressure  variations  are  usually 
associated  with  a  flow  regime  transition.  In  their  work,  the  instabilities 
were  associated  with  transitions  from  bubble-to-plug  flow  and  annu lar- to-s 1 ug 
flow.  In  regard  to  the  general  concept  that  flow  transitions  may  have  an 
important  influence  on  stability,  Isbin(27)^  discussing  a  paper  by  Mendler, 
et  al,  suggested  that,  if  one  provides  step  changes  in  flow  regimes  correspond¬ 
ing  to  specified  conditions,  sustained  oscillations  may  be  generated. 

The  change  in  effective  head  that  occurs  at  the  onset  of  bulk  boiling  can 
be  serious  enough  to  generate  sustained  osc i I lat ions '28  Bergles  and 

Rohsenow^2^j  observed  sustained  oscillations  where  a  transition  takes  place 
between  nucleate  and  Film  boiling.  In  the  latter  case,  film  boiling  begins  in 
a  restricted  area  and  moves  upstream,  followed  by  a  surge  of  water  which  cools 
the  tube.  The  temperature  again  rises  and  the  fluid  flow  decreases  until  a 
transition  is  made  from  nucleate  to  film  boiling  at  a  point  downstream.  In 
tubes  of  small  diameter  (less  than  1/4  in. ID),  the  cycle  is  complete  In 
seconds.  Matzner(30)  also  reports  on  instability  that  could  be  due  to  a 
transition  from  nucleate  to  film  boiling.  In  this  case,  the  instability  took 
the  form  of  a  sudden  transient  pressure  drop  with  an  immediate  recovery. 

While  annular  flow,  mist  flow,  and  fully  developed  slug  flow  are  all 
stable,  both  bubble  flow  and  "developing"  slug  flow  are  transient  regimes  and 
could  be  expected  to  lead  to  flow  instabilities  under  certain  conditions.  An 
instability  in  a  subcooled  water  flow  system  under  conditions  of  zero  gravity 
has  been  reported  by  Pap‘2ll(^*',  which  he  attributes  to  transitions  between 
bubble  and  slug-flow  regimes.  Extremely  severe  flow  perturbations  were 
required,  however,  to  initiate  the  instability.  Under  normal  gravity  conditions, 
the  instability  could  not  be  produced,  even  with  large  perturb.st ions.  This 
has  some  similarity  to  a  phenomenon  reported  by  Radovcich  and  Moissis(^2). 

In  a  no>'m3 1 -gravi  ty  adiabatic  system  of  air  and  water,  they  observed  a  mode 
in  which  Taylor  bubbles  some  50  to  70  tube  diameters  in  length  would  leave  the 
tube  in  v/oat  would  appear  to  have  beer  "developing"  slug  flow.  This  would 
be  fol  lov/ed  by  a  return  to  bubble  flow  and  a  buildup  of  large  bubbles  to 
complete  the  cycle.  In  contrast  to  Papell’s  interpretation,  they  attribute 
this  cycle  to  some  instability  limitation  which  was  not  discussed. 

Vri  view  of  the  earlier  wo''k  of  Isbin,  et  al  ^26,27)^  question  is  not 
definitely  settled,  and  the  possibility  that  this  form  of  instability  may 
appea  ■  cannot  be  definitely  ruled  out. 

Another  unstable  mode  of  operation  that  could  arise  from  a  flow  regime 
transition  involves  the  transition  from  laminar  to  turbulent  flow.  A  similar 
situation  may  be  encountered  in  two-phase  flow  In  the  neighborhood  of  incipient 
boiling  where  bubbles  first  form  on  the  walls.  Such  a  transition  also  could 
take  place  near  burnout  in  the  case  of  highly  subcooled  fluids.  Tripping  of 
the  boundary  layer  from  laminar  to  turbulent  flow  in  the  neighborhood  of 
incipient  boiling  has  been  reported  by  Bergles  and  Rohsenow'^^ ).  Depending 
upon  the  geometry,  this  could  reduce  the  flow  enough  to  cool  the  tube  to  the 
point  where  the  flow  would  become  laminar  and  provide  the  mechanism  for 
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sustained  oscillations.  Desc >•  ipt ions  of  this  mode  of  instability  were  not 
found  in  the  literature.  There  may  be  a  relationship  to  the  flow  oscillations 
in  subcooled  fluids  that  are  associated  with  the  sudden  increase  in  pressure 
drop  taking  place  in  the  neighborhood  of  burnout. 

2.  Negative  Pressure  Drop  vs  Flow  Rate  Derivative 

It  has  been  known  for  a  number  of  years  that,  under  a  wide  range  of 
operating  conditions,  the  pressure  drop  plotted  against  flow  rate  has  a  smooth 
hump  in  an  otherwise  monotonica 1 ly  increasing  curve.  As  illustrated  in 
Figure  then,  there  are  three  possible  flow  rates  for  a  given  pressure  drop. 
If  the  system  should  ope^'ate  at  the  intermediate  point,  it  may  be  seen  from 
the  curve  that  flow  will  increase  with  decreasing  pressure  drop.  This  is  a 
very  unstable  situation,  and  the  system  may  experience  a  very  rapid  increase 
in  flow  rate  until  it  is  operating  at  the  high  flow  rate.  On  the  other  hand, 
the  operating  point  may  shift  rapidly  to  the  low  flow  rate  and  lead  to 
burnout^^^/.  Since  both  high  and  low  flow  rates  are  stable,  the  instability 
consists  of  a  simple  flow  excursion.  This  form  of  instability  has  recently 
been  s'  .^'ied  by  Fraas(*2,33)^  who  calculated  pressure  drop  vs  flow  rates  in 
channe’  ^or  the  purpose  of  establishing  stability  criteria  which,  as  seen  in 
Figure  3  is  assured  if  operating  parameters  are  chosen  to  force  the  system  to 
a  point  where  the  slope  is  definitely  positive.  His  calculations  reveal  that 
the  negative  slope  region  increases  with  the  amount  of  subcooling  and  decreases 
with  the  system  pressure  (see  Figure  4),  confirming  experimentally  observed 
trends.  Of  particular  interest  Is  that  the  Instability  is  impossible  for 
high  pressure  systems  (in  the  neighborhood  of  600  psia  for  water).  Although 
it  is  not  practical  to  have  saturation  conditions  at  the  Inlet,  it  none^  .eless 
is  of  Interest  to  note  that  such  an  operating  point  is  stable  relative  to 
negative  pressure  drops. 

(34 ) 

Recent  experimental  results  of  Collier  have  shown  that,  while 
operation  within  the  negative  slope  region  may  be  impossible,  a  negative  slope 
is  not  a  necessary  condition  of  oscillatory  instability  in  parallel  heated 
channels. 

3.  Compress! bl 1 1 ty 

For  a  number  of  years  authors  have  recogniz€;Ci  that  compressibility  effects 
can  lead  to  both  transient  and  sustained  oscillatiors  in  fluid  flow  systems. 
This  is  quite  evident  for  two-phase  flow  systems  v/here  the  design  results  in 
Impedance  of  the  vapor  phase  exhaust.  It  also  ma/  he  serious  in  single-phase 
flow  systems,  as  was  demonstrated  by  the  theoretical  work  of  Poppendlek  and 
Rickard(35),  in  which  they  showed  that  the  observed  excursions  and  oscillations 
accompanying  reactor  power  surges  could  be  accounted  for  by  considering 
compressibility  and  thermal  expansion.  More  recently,  Lowdermilk,  et  al(23) 
have  called  attention  to  the  fact  that  systems  have  a  strong  tendency  to 
instability  when  designed  in  such  a  way  that  noncondensible  gases  can  become 
entrained  and  trapped.  This  was  demonstrated  by  the  introduction  of  a 
compressible  volume  containing  a  noncondensible  gas  in  the  flow  system  between 
a  throttling  valve  and  -^he  test  section.  Here  the  velocity  of  the  gas  was 
essentially  zero,  so  that  frictional  losses  were  extremely  small  relative  to 
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Figure  3.  Dependence  of  Negative  Pressure 
Drop  Region  on  Subcooling 


feet  of  System  Pressure  Variation  on  the  Negative 
assure  Drop  Regions  at  Constant  Subcooling 
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Figure  4.  Ef 
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those  of  the  vapor  in  two-phase  ^low.  During  burnout  they  observed  flow 
fluctuations  which  Oucurred  at  iower  heat  transfer  rates  as  the  static 
compressible  \olume  was  inc  ased. 

Bergles  and  Rohsenow  29)^  In  3  recent  study  of  burnout  in  tubes  of  small 
diameter,  concluded  that  f low  osc 1 1 lations  caused  by  system  compressibility 
can  greatly  reduce  the  burnout  heat  flux  in  the  subcooled  region,  and  that 
this  instability  is  particularly  difficult  to  avoid  with  sma 1 1 -d i ameter  tubes. 

4.  Parallel  Channel  Coupling 

Because  the  technical  literature  concerning  instabilities  in  parallel 
channel  flow  systems  is  extremely  sparse^^  no  definitive  statements  on  this 
subject  are  possible.  However,  Fraas^*^;  analysis  on  multichannel  stability 
concludes  that  several  possibilities  exist  for  flow  instability  in  parallel 
channels  if  the  individual  channels  operate  in  a  region  that  may  give  rise  to 
flow  instability.  He  notes  that  coupled  oscillations  may  develop  between 
parallel  channels,  inducing  flow  rate  to  increase  In  one  channel  and  diminish 
in  another. 

The  experimental  program  reported  by  Matzner^^^^  is  of  some  interest  in 
connection  with  the  parallel  channel  instability  problem.  Forced  convection 
was  maintained  through  a  bundle  of  19  parallel  electrical  heaters  wh’  1  were 
packed  close  enough  to  very  roughly  approximate  a  number  of  parallel  lannels. 
At  high  pressures,  only  slight  pressure  fluctuations  were  observed.  On 
approach  to  burnout,  pressure  transients  were  attributable  to  pc-sible  flow 
regime  transitions  as  discussed  above.  The  inside  rod  did  not  t  dergo  any 
temperature  instability,  while  the  adjacent  six  inner  rods  showed  worse  per¬ 
formance  than  the  outer  12  rods.  Of  particular  interest  was  the  fact  that 
the  pressure  transient  took  the  form  of  a  pressure  drop  that  did  not  exceed  4 
to  12  psi  at  an  operating  pressure  of  1000  psi.  This  was  followed  by  an 
Immediate  recovery  to  slightly  over  the  initial  pressure. 

Pi scuss ion 

The  literature  review  summarized  above  revealed  far  more  problems  than 
it  solved  and  led  to  the  obvious  conclusion  that  the  two-phase  flow  stabi  lity 
problem  is  extremely  complex.  There  are  many  different  types  of  instabilities; 
each  may  exhibit  more  than  one  behavior  and  is  affected  by  numerous  parameters, 
including  fluid  properties,  geometry,  and  imposed  boundary  conditions. 
Instabilities  may  exist  with  steady  imposed  conditions,  or  may  be  caused  by  a 
change  in  the  nominal  operating  condition  of  the  heat  exchanger. 

It  does  not  appear  possible  to  derive  one  equation  covering  all  possible 
instabilities;  each  type  of  instability  and  set  of  boundaiy  conditions 
generally  must  be  analyzed  separately.  iVhiia  most  types  have  been  observed 
empirically,  usually  by  unfortunate  accident,  few  have  been  correlated 
analytical ly.  Most  of  the  analyses,  while  predicting  qualitatively  the  effect 
of  some  of  the  variables,  have  been  unsuccessful  in  predicting  the  quantitative 
results  which  were  experimentally  observed.  In  addition,  none  of  the  analyses 
ha«,  been  verified  thoroughly  enough  by  experiment  to  warrant  its  application 
without  subsequent  experimental  confirmation. 
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It  is  fel^  that,  in  a  once- through,  high  exit  quality  vaporizer,  thp  n’os  t 
prevalent  class  of  instability  is  that  resulting  from  negative  derivative  of 
pressure  drop  with  respect  to  flow  rate.  The  analysis  of  this  instability 
by  Fraas  (33)  has  much  to  recommend  it  since  it  is  simple,  reveals  the  phys¬ 
ical  importance  of  the  variables,  and  has  predicted  the  effect  of  pressure. 
Subcooling,  and  heat  flux  on  flow  stability  which  has  been  observed  in  practice. 

In  summary,  the  state  of  the  art  is  not  far  advanced;  there  is  lack  of 
agreement  about  the  effect  of  many  variables;  most  conclusions  about  two- 
phfise  flow  stability  are  qualitative  (e.g.,  high  pressure  is  good)  and  there 
is  no  general  analytical  approach  which  shows  promise  of  describing  flow 
stabi 1 i ty  behavior, 

SIMULATION  CRITERIA 

Because  of  the  complexity  of  the  two-phase  flow  stability  prob’em  and 
the  lack  of  an  available,  dependable,  general  analysis  of  the  problem,  the 
decision  was  made  to  resolve  the  SNAP-50/SPUR  boiler  flow  stability  problem 
experimentally.  A  through  investigation  of  tie  effect  of  all  the  important 
parameters,  using  potassium  as  the  test  fluid,  was  not  economically  feasible. 
Also,  it  was  judged  very  desirable  to  use  a  fluid  which  would  allow  the  con¬ 
struction  of  a  transparent  boiler  so  that  visual  observations  would  be  pos¬ 
sible;  potassium  does  not  satisfy  this  requirement.  It  was  therefore  neces¬ 
sary  to  develop  simulation  criteria  which  would  provide  a  rational  basis  for:, 
selecting  a  test  fluid,  designing  a  test  loop,  and  planning  a  test  program 
whose  results  would  be  applicable  to  the  SNAP-50/SPUR  potassium  boiler. 

The  literature  rev'ew  revealed  a  large  number  of  parameters  affecting 
the  two-phase  flow  stability  problem.  Those  discussed  by  various  investi¬ 
gators  are 
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The  expense  of  experimentally  invest igar i nn  a  range  of  — — 
of  the  above  parameters  is  prohibitive,  even  with  a  conventional  fluid.  Be¬ 
cause  of  this, and  the  necessity  for  applying  the  experimental  results  to 
potassium,  it  was  judged  necessary  to  reduce  the  important  physical  variables 
to  a  complete  set  of  independent,  dimensionless,  similarity  parameters. 


The  application  of  dimensional  analysis  to  heat  transfer  and  fluid  me¬ 
chanics  problems  is  very  common.  It  serves  two  very  important  functions  in 
experimental  investigations:  (l)  the  number  of  independent  variables  is 
reduced  because  there  are  fewer  independent  dimensionless  groups  than  physi¬ 
cal  variables,  and  (2)  simulation  of  the  actual  problem  is  possible  with  more 
convenient  geometry  and  test  fluids.  The  theory  and  methods  of  dimensional 
analysis  are  discussed  by  McAdams  (39)  and  other  authors  of  heat  transfer 
textbooks.  The  method  used  here  to  define  dimensionless  groups  was  to  com¬ 
bine  (l)  nondimens ional iz i ng  of  the  controlling  fundamental  equations,  (2) 
past  experience  with  the  same  physical  variables,  (3)  and  intuition.  The 
nondimensional izing  procedure  was  simplified  because  Fraas  (12)  presented 
his  results  in  dimensionless  form. 


In  calculating  the  pressure  drop-vs-flow  rate  relations,  it  becomes 
apparent  that  the  importance  of  subcooling  is  determined  by  the  ratio  of  (l) 
the  heat  added  to  raise  the  fluid  to  the  saturation  temperature  to  (2)  the 
heat  rer;uired  for  vaporization.  Therefore,  the  dimensionless  parameter  that 

measures  the  effect  of  subcooling  is  C  AT/hr  •  This  same  parameter  appears 

P  ^9 

if  the  energy  balance  for  a  boiler  is  expressed  in  dimensionless  form.  It 
also  has  been  revealed  in  analysis  of  boiling  heat  transfer;  and  in  soma 
places  is  called  the  Boiling  Number. 


As  discussed  above  in  the  literature  review,  pressure  level  is  Important 
primarily  because  of  its  effect  on  density  ratio.  For  a  constant  friction 
factor,  the  ratio  of  liquid  to  vapor  pressure  drop  is  proportional  to  the 
density  ratio,  which  is  therefore  a  measure  of  potential  pressure  drop  fluc¬ 
tuation.  For  these  reasons,  and  because  previous  investigators  have  used 
1 iquid-to-vapor  density  ratio  as  an  independent  variable,  it  was  selected  as 
one  of  the  dimensionless  similarity  parameters. 


Pressure  drop,  flow  rate,  inlet  velocity,  and  mass  velocity  were  chosen 
separately  for  emphasis  by  various  investigators.  However,  they  are  all  a 
function  of  the  dimensionless  pressure-drop  ratio,  (pressure  drop  divided 
by  pressure  level,)  which  is  proportional  to  Mach  number  and  L/D.  The  system 
geometry  Is  defined  by  the  length-to-diameter  ratio. 

The  heat  balance  on  a  boiler,  expressed  in  dimensionless  form,  shows 

(J/A 

that  the  dimensionless  heat  flux,  ^iT"  ^  ®  function  of  the  dimensionless 

fg 

subcooling  parameter,  exit  quality,  and  L/0.  Therefore,  exit  quality  was 
chosen  as  the  final  independent  variable. 


Subcooled  liquid  pressure  drop,  expressed  in  dimensionless  form  by 
dividing  it  by  the  two-phase  pressure  drop,  was  selected  as  the  dependent 
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d  ('-.lens  ion  less  parameter.  The  results  of  Fraas  Lowdermilk  'JO*,  and 

others  showed  that  tv/o-phase  flow  always  could  be  stabilized  with  a  great 
(’noug.h  pressure  drop  in  the  subcooled  liquid  region.  As  a  result,  the  de¬ 
cision  was  made  to  solve  the  two-phase  flow  stability  problem  in  the  SNAP- 
50/SPUR  boiler  by  itsertmg  orifices  at  the  entrance  of  each  bo‘ ler  tube. 

The  lest  program  described  in  the  following  section  was  designed  to  provide 
spic'fic  information  pertaining  to  the  magnitude  of  the  pressure  drop  in  the 
subcooled  liquid  region  required  to  provide  stable  flow. 

In  Summary,  the  physical  variables  were  reduced  to  a  complete  set  ot 
s '  j<  independent,  dimensionless  parameters,  which  can  be  expressed  as  fol'o.vs 


AP 


The  purpose  of  the  test  program  was  to  determine 
liquid  pressure  drop  ratio  required  to  produce  stable 
five  independent  parameters  on  the  right-hand  side  of 
a  range  abound  SNAP-50/SPUR  boiler  test  and  operating 


experimentally  the 
flow  for  values  of  the 
the  above  equat'on  cr//er’n{ 
COnd i t i Ons 
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SECTION  -3 


TEST  PROGRAM 


OBJECTIVE 

The  primary  objective  of  this  program  was  to  obtai^^  Sufficient  infor¬ 
mation  to  design  the  SNAP  50/SPUR  boiler  so  that  stable  two-phase  flow  will 
occur.  The  boiler  is  a  tube-and-shei 1  up't  with  twisted  tapes  inserted  in 
the  tubes  to  produce  swirl  which  wi 1 i  centrifuge  entrained  droplets  to  the 
tube  wall.  The  design  approach  was  to  place  an  orifice  'n  the  subcooled 
liquid  region  at  the  entrance  to  each  boiler  tube.  The  liquid  pressure  drop 
across  the  orifice  must  be  made  high  enough  to  dominate  the  possible  vari¬ 
ations  in  two-phase  flow  race  that  are  possible  with  a  constant  header-to- 
header  pressure  drop.  The  ratio  of  orifice  pressure  drop  to  boiling  two- 
phase  pressure  drop  is  considered  to  be  the  controlling  d'mens ionless  param¬ 
eter  and  must  be  maintained  above  a  critical  va'ue  to  provide  stable  flow. 
Since  the  problem  was  too  complex  to  analyze,  this  test  program  was  developed 
to  determine  the  pressure-drop  rat  o  reqi,jred. 

A  major  objective  for  later  work  w'th  the  test  loop  des', 'ibed  below  is 
a  study  of  the  condensing  process  in  parallel  tubes.  Since  the  same  loop, 
type  of  instrumentation,  and  observations  will  be  util'zed  in  the  condensing 
studies  as  in  the  boiling  studies,  the  objectives  and  requirements  of  this 
later  work  necessarily  affected  loop  design. 

Due  to  the  complexity  of  the  two-pnase  flow  process,  more  than  one  plaus¬ 
ible,  but  contradictory,  explanation  often  can  be  found  to  rationalize  un¬ 
expected  experimental  results,  for  thi?  reason,  it  is  essential  that  research 
on  iwo“ph3S('.-  flow  include  visual  observations  whenever  possible.  The  follow¬ 
ing  partia-  list  gives  examples  of  the  type  of  information  that  only  visual 
observations  can  provide. 

^  Determine  the  flow  patterns  that  occur  during  forced-convection 
vaporization,  wth  and  without  tube  inserts. 

2.  Determine  if  a  substantial  amount  of  liqu'd  flows  along  the  twisted 
tape  after  the  wall  is  dry,  and  how  the  liquid  breaks  off  the  end 
ot  the  tape. 

3.  Oete  rrri  ne  the  length  beyond  the  ena  of  a  twisted  tape  that  the  swirl 
vs^vch  was  induced  by  the  tape  persists. 

4.  Establish  the  mechanisms  for  the  various  flow  regime  transitions. 

5.  Establish  the  flow/  patterns  and  flow  maldistribution  that  occurs 
in  parallel  condensing  tubes 

6.  Study  the  fluid  dynamics  at  the  condenser  liquid-vapor  interface. 

7.  Establish  the  visual  characterist ics  of  boiler  and  condenser  flow 
I  ns  tab i 1 i t ies . 
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The  test  program  was  designed  to  achieve  the  general  objectives  described 
above  for  values  of  the  controlling  dimensionless  parameters  developed  in  the 
preceding  section,  including  the  range  of  importance  to  the  SNAP  50/SPUR 
boiler.  Table  I  presents  the  current  values  of  the  five  dimensionless  simi¬ 
larity  parameters  in  the  design-point  boiler  and  in  a  19-tube  prototype  test 
boiler. 


TABLE  1 

VALUE  OF  DIMENSIONLESS  SIMILARITY  PARAMETERS  IN  SNAP  50/SPUR  BOILER 


Design  Point  Boi ler 

19-Tube  Test- 

c  AT/h, 

P  ^9 

0.18 

Variable 

200 

400 

£.P,p/P 

0.  1 

0.4 

Exit  Quality,  percent 

100 

Variable 

L/D 

150 

200 

In  summary,  the  test  program  and  test  loop  were  designed  to  determine 
the  pressure  drop  in  the  subcooled  liquid  region  of  the  boiler  required  to 
provide  stable  two-phase  flow  for  a  range  of  conditions  including  those  pre¬ 
sented  in  Table  1;  it  was  necessary  that  this  be  done  with  a  transparent  test 
loop. 

FLUID  SELECTION 

A  survey  of  the  compatibility  of  the  characteristics  of  common  fluids 
with  the  test  objectives  of  this  program  led  to  selection  of  Freon  113  as 
the  test-loop  working  fluid.  This  inexpensive,  easily  available  fluid  is 
relatively  nonhazardous  and  has  well  established  physical  properties. 

To  decrease  the  hazards  involved  in  visually  observing  the  two-phase 
flow  processes  at  close  range,  it  was  desirable  to  run  the  tests  near  atmos¬ 
pheric  pressure  and  temperature.  The  saturation  temperatures  of  Freon  113 
are  80®Fand  120®F,  respectively  for  1 iqu id-to- vapor  density  ratios  of  400  and 
200.  The  saturation  pressures  at  these  conditions  are  7  and  15  psia,  as  shown 
in  Figure  5.  While  both  saturation  pressure  and  temperature  are  somewhat 
lower  than  the  most  desirable  values,  they  are  compatible  with  the  antici¬ 
pated  test  program  and  the  requirement  for  maintaining  the  boiling  temperature 
sufficiently  above  that  of  common  cooling  fluids  to  ease  the  design  of  the 
condensing  side  of  the  test  loop. 

Another  advantage  of  Freons  in  general,  and  Freon  113  in  particular,  is 
the  relatively  low  value  of  60  to  70  Btu  per  lb  for  the  heat  of  vaporization. 
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Figure  5  Freon  113  Properties 
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dimensionless  sobcooling  paramete-',  s-nce  only  small  cnanges  -n  entrance 
temperature  are  'equired.  In  addition,  both  the  heat  and  the  total 

quantity  of  heat  -eqo-red  to  vaporize  the  flu-d  completely  are  proportional 
to  the  heat  of  vapor ■ zat ion .  low  values  for  these  parameters  ease  the 
test  loop  design  prop 'em 

^  finally,  freon  systems  are  common  at  AiResearch;  therefore,  a  great  deal 
o.  e-xperience  existed  for  the  design  and  proper  instrumentation  of  a  reliable, 
leakproof  freon  test  loop. 

RUN  SPECIf ICATION 

Of  the  five  parameters  given  ir>  fable  i,  't  was  necessary  to  test  a 
range  of  values  of  only  the  first  three,  s ■ nee  the  ex , t  qua'-ty  was  kept  near 

percent  in  all  r,ns,  and  the  test  boiler  was  designed  with  a  constant  L/D 
rat  I o  of  1 60. 

fhe  dimensionless  Subcooling  parameter  has  the  required  value  of  0  IS 
with  Freon  113  when  the  entrance  subcooimg  is  SO'-f.  To  determine  the  effect 
of  variations  in  subcooling,  values  25°f,  50°f,  and  l00°f  were  tested. 

corresponding  to  density  ratios  of 
00,  200,  and  50,  respectively  were  tested.  At  each  of  these  pressures, 

P  essure  drops  o^  1,  2,  and  3  psi  were  run,  y^eld-ng  a  normalized  pressure 
roD  range  of  0.05  to  0-3  The  lowest  pressure  of  10  psia  .s  higher  than 
the  value  of  7  ps.a,  whicn  corresponds  to  the  l9-iuDe  boiler  dens-ty  ratio, 
owever,  It  was  the  minimum  pressure  that  allowed  substantial  subcoolinq  to 
e  tested,  m  the  loop  desenped  be'ow,  due  to  the  difficulty  of  coolmq  the 
oi  er  in  et  to  much  below  50  F.  the  pressure  range  provides  a  2  to  I  vari- 
boiler'"  pressure  simulates  me  design-point 


The  individual  test  conditions  r^n  for  each  type  of  tube  nsert  as 
isted  in  Table  2,  include  most  combinations  of  the  following  properties 

^nur  ■  20  ps-a 


T  .  T 
sat  in 


-  100,  50,  and  25'''f 

=  1,2,  and  3  psia 
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TABLE  2 


BOILER  TEST  PROGRAM 


RUN 

P  , 
out 

(psia) 

iPip 

(psia) 

sat  ID 
(°F) 

^out 

(°F) 

T. 

1  n 

/Op) 

1 

20 

1 

100 

135 

40 

2 

20 

1 

50 

135 

90 

3 

20 

1 

25 

135 

l!5  1 

4 

20 

2 

100 

135 

40  1 

5 

20 

2 

50 

135 

90  j 

6 

20 

2 

25 

135 

!  15  i 
) 

7 

20 

3 

100 

135 

43  i 

J 

8 

20 

3 

50 

135 

95  i 

f 

9 

20 

3 

25 

135 

120  ! 

!  1 

15 

1 

100 

120 

20  1 

!2 

15 

1 

50 

120 

70  1 

13 

15 

1 

25 

120 

95  ; 

14 

15 

2 

100 

120 

25 

15 

15 

2 

50 

120 

75 

16 

15 

2 

25 

120 

!00  i 

17 

15 

3 

100 

!20 

30 

18 

15 

3 

50 

120 

30 

19 

15 

3 

25 

120 

105 

24 

10 

1 

25 

100 

75 

25 

!0 

1 

50 

too 

5u 

26 

10 

2 

25 

100 

80 

27 

10 

2 

50 

100 

55 

28 

10 

3 

25 

too 

35  ■; 

29 

_ 

10 

3 

50 

loo 

60 

23 


SECTION  4 


APPAJ^ATUS 


TEST  LOOP 

A  sche'ifitic  of  the  Preon  f  low-stab i  i  <  ty  test  loop  -s  shown  in  Pigu'^e  6. 
figures  1,  (I,  and  9  a'^e  photog'-aphs  of  the  loop,  showmg  the  DO' ler^  test 
condenser,  au/iiliary  equipment,  and  •  nsfumentat  ion  pane*  The  two  compon¬ 
ents  stuciied  were  boile'"  and  test  condenser,  both  of  these  were  made  from 
transpa'^ent  materials,  as  was  much  of  the  rest  of  the  loop  As  mentioned 
previously,  the  primary  parameter  for  this  study  was  the  • i qu i d-to- two-phase 
pressure-drop  ratio. 

A  va’iabie  orjPice  (Hoke  280  sertes  metering  valve),  as  shown  in  Pigure  9 
was  installed  near  the  entrance  of  each  boiler  tube  The  Preon,  after  leaving 
the  boiler,  passed  through  a  glass  water-cooled  heat  exchanger.  The  main 
use  of  this  heat  exchanger  was  'n  conjunction  with  condensing  tests,  the  qual¬ 
ity  to  the  condenser  could  be  adjusted  at  the  cooler  wh'le  maintaining  con¬ 
stant  boiler  conditions. 

The  mam  condenser  was  capable  of  condensing  aM  the  "^eon  vaporized  in 
the  boiler.  The  cooling  medium  was  air,  passing  counterflow  to  the  Preon. 
Liquid  coolants  also  may  be  used. 

A  versatile  mounting  for  the  high-speed  camera  was  fabricated  above  the 
boiler  as  shown  in  Pigure  9.  The  mounting  was  such  that  the  camera  could 
film  any  tube  and  any  pos'tion  along  the  tube 

Liquid  Preon,  leaving  the  condensmg  section  of  the  loop,  passed  through 
a  subcooler  (chilled  by  water  or  tr < ch loroethy lene)  before  entering  the  pump. 
The  pump  was  an  exblos 'On-proof,  seal-less,  centrifugal  Series  G  Chempump. 

A  wraparound  water-cooling  jacket  was  incorporated  for  pump  cool  mg 

A  bypass  loop  and  control  valve  were  incorporated  around  the  pump.  Ail 
tests  were  conducted  with  the  bypass  control  valve  m  a  cmsed  position.  A 
Greer  bladder  accumu lator,  pressur i zed  with  wate",  was  ut'l'zed  for  system 
makeup  and  pressure  level  control.  Host  of  the  tests  were  conducted  at  pres¬ 
sure  levels  that  allow  the  solenoid  vaive  to  be  'n  a  dosed  position. 

COMPONENTS  AND  INSTRUMENTATION 

The  boiler  which  'S  shown  in  Pigu-e  10,  consists  of  five  quartz  tubes 
(6-mm  ID  ,  8-mm  Ol),  and  53-in.  length)-  A  0.63  m  by  3.57  m  'ectangular 
cross-section  pyrex  duct,  42-1/2  in  m  length,  was  bu  ■  1 1  to  Surrour.d  the 
tubes  so  that  hot  am  would  flow  outside  the  tubes  and  counterflow  to  the 
Freon.  The  air  flow  rate  was  limited  to  18  lb  per  mm,  with  a  maximum  inlet 
temperature  of  800°F.  Figure  I  I  >5  a  schematic  of  the  boiler  showing  those 
quantities  that  were  measured. 
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Figure  7.  Freon  Flow  Stability  Loop 
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The  test  condenser  consisted  of  30  quartz  tubes  the  same  dimensions  as 
the  bo  i  lei"  tubes.  .A  rectangular  cross-sect  ion  FyreA  glass  duct,  42-i/2  in. 
in  length  was  constructed  around  the  tubes  to  serve  as  the  coolant  flow  pas¬ 
sage. 

An  electric  heater  (6-kw  capacity)  was  provided  to  control  the  Freon 
inlet  temperature  to  tne  boiler.  An  electrically  operated  Series  715  Flow- 
Ball  Valve  was  used  to  control  the  liquid  flow  rate  to  the  boiler,  which  was 
measured  with  a  turbine-type  Cox  flowmeter. 

Ashcroft  Laboratory  test  gauges  were  used  to  measure  pressures  in  the 
Freon  loop.  Transducers  (Pace  Engineering  Model  PID  and  Model  PIG,  respec¬ 
tively)  were  used  to  i-ecord  the  liquid  (orifice)  pressure  drop  and  the  total 
boiler  pressure  drop.  A  standard  calibrated  orifice  was  used  for  air-flow 
measurements.  Inlet  and  outlet  temperatures  on  all  heat  exchangers  were  mea¬ 
sured  with  i ron-constantan  thermocouples  and  recorded  on  a  Brown  48-pO'nt 
temperature  indicator.  Toggle-action  Dahl  valves  were  used  throughout  the 
loop  to  control  rhe  Freon  flov/  rate  to  the  various  heat  exchangers.  Figures 
12  and  13  are  photographs  of  the  instrumentation  and  control  panel. 

Table  3  lists  components  used  'n  the  loop.  Glass  piping  and  fittings 
were  used  when  convenient. 

PHOTOGRAPHIC  EQUIPMENT 

A  WF-3  ^astax  camera,  mounted  horizontally  above  the  test  section  as 
shown  in  Figure  14,  obtained  vertical  shots  by  use  of  a  2  by  2  prism  Single¬ 
tube  c  :>seups  were  shot  with  a  152-mm  lens  in  the  camera  and  a  supplementary 
6-in.  Mero  Ektar  lens  mounted  horizontally  below  the  prism. 

Shots  of  all  five  tubes  were  taken  with  a  50-mm  lens  in  the  Fastax.  Ad¬ 
ditional  five-tube  shots  were  obtained  by  using  a  Cine  Kodak  Special  camera 
with  either  a  25-mm  or  50-mm  lens. 

Ektachrome  ERB  (ASA  125)  fil  m  was  used  in  the  Fastax  camera  and  Ektachrome 
Commercial  Type  72S5  (ASA  25)  film  was  used  in  the  Kodak  camera. 

Lighting  was  provided  by  two  Colortran  Quartz  i i te  1000-w,  1 lO-v  lamps 
located  below  the  boiler,  which  gave  50,000  to  100,000  ft-candles  of  light. 

The  cameras  were  stopped  down  to  provide  the  maximum  depth  of  field  possible 
with  the  above  lighting.  When  the  Fastax  camera  was  operating  at  its  maxi¬ 
mum  film  speed  of  BOOO  frames  per  sec  in  the  single-tube  shots,  the  resulting 
depth  of  field  was  approximately  0  10  in. 
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TABLE  3 


FREON  LOOP  COMPONENTS 


Item  Number 
on  Fiqure  6 


Cotuponent 


Description 


Boi ler 


5  quartz  tubes,  53-in,  length,  6-mm  ID, 
8-mm  OD;  Pyrex  duct  built  over  tubes 
for  a  i  ■"  flow 


Test  Condenser 


Cooler 


Main  Condenser 


Sub-Coo ler 


Accumu la  tor 


Solenoid  Valve 


Pump 


Electric  Heater 


Liquid  Control  Valve 


Variable  Orif'ces 


Water  Flow  Meters 


Freon  Control  Valves 


Pressure  Gauge 


30  quartz  tubes,  53-in.  length,  6-rnm  ID, 
B-rnm  OD;  Pyrex  duct  built  over  tubes 
for  air  flow 

QVF  2/3i  gl  ass  condenser 

Q,VF  6/15  glass  condenser 

Q.VF  2.3^  glass  condenser 

lO-ga!  Greer  transfer  barrier  accumu¬ 
lator  with  a  neoprene  bladder,  coated 
shell  (phenolic  resin  sprayed  and 
baked),  300-series  stainless  steel 
fittings. 

Sporlan  solenoid  valve,  l-in.  FPT, 
direct-acting,  24-s/  d-c  coil 

Series  G  Chempump,  Model  GA  I  1/2 

Thermolab  duct  assembly  preheater, 

1  I 5-v,  OO-cyc !e 

Series  716  flow- ball  vaive,  electricaMv 
operated. 

Hoke  2 JO  series  metering  valve,  globe 
pattern,  3I6SS  valve,  l/8-in.  orifice 

Sho  rate  ’’50"'  flowmeter,  Model  1350, 
size  8,  chrome  plated,  brass  const''ucted 

Dahl  valves;  panel  mounting;  toggie 
action;  Teflon  seat 

Ashcroft  Laboratory  test  gauge,  60  ps i 
to  30- i n .  vacuum 


Temperature  Indicator 


Brown  temperature  indicator,  48-point, 
IC  thermocouple,  i I 5-v  60-cycle,  I00°F 
to  i00C°F  scale.  M i nneapo 1 i s'Honeywe 1 1 


SECTION  5 


TEST  PROCEDURE 


RUN  CRITERIA 

The  tolerances  allowed  on  the  nor.inal  test  conditions  for  each  run  pre¬ 
sented  in  Table  2  were  ±0.5  psia  on  I’  ,  ±5®F  on  T.  ,  and  ±0.2  psi  on  APto* 

oil'"  in'  TP 

Each  run  (set  of  test  conditions)  conn  sted  of  obtaining  data  at  nominal  valve 

settings  of  2, 3, 4, 5, 8,  and  15  (pressu  *<>  drop  across  the  valve  is  a  function 

of  the  valve  setting,  with  0  the  full  closed  posit icn).  The  i uns  were  started 

with  the  valve  setting  of  two,  which  always  gave  stable  flow.  The  air  and 

Freon  flow  rates  were  adjusted  to  /ie  <  the  soecifled  boiling  pressure  drop 

(P^  ,  -  P  and  an  exit  quality  cltue  to  1.0. 

tube  out 

The  valve  settings  on  the  fiva  tutes  were  adjusted  so  that  the  average 
value  of  the  tube  pressure  read  (rr or  eich  gauge  was  the  same  within  ±0.2  psi. 
This  provided  uniform  flow  in  each  tuba-  Other  checks  on  the  uniformity  of 
conditions  in  the  five  tubes  were  the  visual  appearance  of  the  flow  ir.  the 
Freon  tubes  on  either  end  of  t**?  bo’Uir,  and  the  uniformity  of  the  transducer 
traces  from  the  five  tubes. 

Air  at  temperatures  up  to  (')0®F  w*  >  used  to  vaporize  the  Freon.  The  air 
flov/  rate  and  temperature  were  adjus’,e<'  so  that  the  temperature  of  the  Freon 
at  the  boiler  exit  was  approximately  3(  °F  higher  than  the  saturation  temper¬ 
ature  corresponding  to  the  boiler  exit  :rcssure.  The  boiler  was  operated  at 
the  highest  exit  quality  possible  wlhi'  the  preceding  restriction  on  boiler 
exit  temperature. 

The  data  points  at  valve  .sttiniK  greater  than  two  were  taken  after 
changing  only  the  valve  settings,  fills  sometimes  tended  to  cause  changes  in 
Freon  exit  quality;  this  was  contiij  lad  by  adjusting  the  air  temperature. 
Tolerances  on  the  primary  variables  For  any  set  of  six  points  were  as  follows; 

Freon  flowmeter,  ^4  cps 

Freon  exit  pressure,  ±0. j  ps i 

Freon  tube  pressure,  ±0.1  p  < 

Freon  inlet  temperature.  y3®F 

The  boiler  pressure  drop  was  within  ±0.2  psi  of  the  specified  value  and 
was  held  as  constant  as  possible  during  each  set  of  six  points,  which  make 
up  a  run. 

A  t/pical  set  of  data, which  corresponds  to  Run  26,  is  presented  in 
Figure  15.  The  data  processing  procedure  described  below  uses  these  data  as 
the  example. 
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HEAT  LOSS  CALIBRATION 


A  heat  loss  calibration  was  performed  initned i ate  1  y  before  and  after  each 
scries  of  runs  by  obtaining  the  following  data  with  no  Freon  in  the  boiler; 
for  air  inlet  temperatures  of  BOO^F,  600°F,  and  400°F,  the  air  exit  temper¬ 
ature  was  measured  for  air  flow  rates  of  9,  12,  15,  and  18  lb  per  min.  These 
data  were  taken  very  carefully,  allowing  sufficient  time  (15  to  30  min)  for 
the  exit  temperature  to  reach  steady  state  after  each  change  in  conditions. 

The  results  of  the  heat  loss  calibration  indicated  that  a  significant 
part  of  the  energy  loss  was  due  to  air  leakage.  An  energy  balance  on  the 
boiler,  assuming  that  the  air  leakage  was  proportional  to  the  air  pressure, 
revealed  that  the  energy  loss  due  to  leakage  was  proportional  to  t^’  i  'oduct 
of  the  air  pressure  and  the  air  temperature  chanae  in  the  boiler.  ini^  loss 
was  added  to  the  heat  leak  itself,  which  was  proportional  to  the  temperature 
difference  between  the  hot  air  and  the  atmosphere.  An  equation  for  the  energy 
loss,  consisting  of  the  two  terms  defined  above,  was  developed  from  the  heat 
loss  calibration  data  which  correlated  all  the  results  within  ±10  percent. 

This  equation  remained  valid  even  after  the  boiler  was  taken  apart  and  re¬ 
assembled  during  the  test  program. 
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SECTIOM  6 


DATA  PROCESSING 


The  data  obtained  for  each  run  included  the  information  recorded  on  the 
data  sheet  for  the  six  valve  settings  and  the  corresponding  transducer  strip 
charts.  A  typical  complete  set  of  data  for  Run  26  with  full-length  twisted 
tape  inserts  is  presented  in  Figures  15  to  21  and  will  be  used  to  illustrate 
the  data  processing  procedure.  The  transducer  strip  chart  calibration  infor¬ 
mation  is  presented  in  Table  4. 

The  data  from  the  data  sheet  shown  in  Figure  15  were  entered  on  the  data 
processing  sheet  as  shown  in  Figure  22. 


TABLE  4 

TRANSDUCER  STRIP  CHART  CALIBRATION 


Trace 
■eak  No. 

Measurement 

Range 
(psi  ) 

Reference 

Point 

(in.) 

Cal i brat  ion 
V  p  .>  i  /  i  n .  ) 

3 

Reference 

4 

iP^,  tube  . 

-2  to  +.0 

1.5 

8 

5 

iP. ,  tube  2 

-2  to  +10 

2.0 

8 

6 

tube  3 

-2  to  -^10 

2.5 

8 

7 

^P^,  tube  4 

-2  to  -10 

3.0 

8 

e 

AP^,  tube 

-2  to  +10 

3.5 

8 

9 

Boiler  P. 

b 

-2  to  +10 

4.0 

4 

10 

Condenser  P 

-2  to  -!0 

4.5 

4 

16 

BrMer  inlet  P. 

1  n 

0  to  50 

6.5 

10 

iS 

Reference 

6.5 

The  exit  quality  was  calculated  on  a  computer.  The  net  heat  input  was 
obtained  by  subtracting  the  calibrated  energy  loss  from  the  energy  decrease 
of  the  air  is  it  flowed  through  the  boiler.  The  Freon  exit  quality  was  then 
calculated  from  an  energy  balance  using  the  net  heat  input. 
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igure  16  Transducer  Strip  Chart  Data,  Pull  Tape,  No.  26,  Condition  No. 


Bo i ler 


Cure  *6  Strip  Chart  Data. 


Boiler  4.? 


igure  17. 


Condenser 


igure  20.  Transducer  Strip  Chart  Data.  Tull  Tape,  Run  No  26,  Condition  No. 


5-TiJBE  FREOH  BOILER  FLOW  STABILIl  Y  DAlA 


RUN  NO. _ 26 _  Date _ 4/30/64 _ 

Tube _ Full-Length  Twisted  Tape  Log  Sheet  No. _ 7^ 
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Figure  22.  Data  Processing  Sheet  for  5-Tube 
Boiler  Flow  Stability  Data 


46 


te*?i 


The  frequency  and  amplitude  of  the  flow  oscillations  were  measured  on 
the  transducer  strip  charts  and  entered  at  the  bottom  of  Figure  22;  the 
overall  boiler  pressure  drop  trace,  break  No.  9,  was  used  for  the  amplitude 
measurement.  The  amplitude  in  psl  was  obtained  by  multiplying  the  measured 
amplitude  by  the  A>psi>per«in.  scale  factor.  After  the  amplitude  was  nor¬ 
malized  by  dividing  it  by  the  boiler  two-phase  pressure  drop,  the  results  were 
plotted  as  shown  in  Figure  23.  The  critical  liquid  pressure  drop  Is  located 
by  the  knee  of  the  curve. 

The  dividing  line  bet%^en  stable  and  unstable  conditions  was  much  sharper 
than  expected.  It  was  actually  possible  to  pick  out  the  critical  liquid 
pressure  drop  by  simply  surveying  the  transducer  strip  charts.  In  the  example 
given,  the  flow  clearly  becomes  unstable  between  conditions  2  and  3.  This 
is  confirmed  by  the  uncertainty  in  the  tube  pressure  measurements  recorded 
by  the  technician,  which  exhibit  a  factor-of-ten  increase  between  these  two 
points.  It  was  possible,  therefore,  to  pick  out  the  critical  liquid  pressure 
drop  directly,  without  going  through  the  steps  represented  by  Figures  22  and  23. 
In  this  case,  the  liquid  pressure  drop  of  condition  two  was  selected  as  that 
required  to  provide  stable  flow  for  the  conditions  of  run  26. 

The  oscillation  amplitude  of  the  overall  boiler  pressure  drop  was  about 
half  the  amplitude  in  the  individual  tubes.  Indicating  a  damping  effect  in 
the  manifolds  where  the  oscillations  of  the  individual  tubes  were  superimposed 
upon  one  another.  The  oscillation  frequency  was  greatest  when  the  valves 
were  wide  open  at  a  nominal  setting  of  15.  It  is  this  maximum  frequency  that 
is  reported  for  each  run  in  the  following  section. 
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RESULTS 


The  experimental  results  are  presented  on  the  following  pages  in 
Tables  5  and  6.  Table  5  includes  the  results  with  twisted  tape  inserts 
running  the  entire  length  of  the  tubes.  The  results  presented  in  Table  6  were 
obtained  with  twisted  tape  inserts  which  extended  only  three-quarters  of  the 
boiler  heated  length;  they  stopped  approximately  iO  in.  from  the  end  of  the 
heated  length.  A  tape-twist  ratio  of  3  (480-degree  twist  in  a  length  equal 
to  three  tube  diameters)  was  used  in  all  tests.  The  tapes  were  relatively 
loose  fitting  and  could  be  pulled  cr  pushed  through  the  tubes  by  hand. 

No  plain-tube  results  are  presented  because  these  earlier  tests  were 
run  with  valve  settings  of  2,  5,  B,  10,  and  15.  This  did  not  provide  close 
enough  definition  of  the  critical  liquid  pressure  drop  since,  in  most  cases, 
a  valve  setting  of  3  or  4  corresponded  to  the  critical  condition.  If  the 
plain  tube  tests  are  completed  by  obtaining  results  at  these  valve  settings, 
the  results  will  be  reported  separately. 

The  results  presented  in  the  tables  include  those  measurements  required 
to  define  the  state  of  the  two  fluids  at  the  highest  valve  setting  which 
gave  stable  flow.  The.  transducer  measurements  corresponding  to  a  valve  setting 
of  15  also  are  presented.  The  amplitude  given  is  that  for  the  overall  boiler 
pressure  drop,  which  generally  is  half  the  amplitude  of  the  individual  tube 
oscillations.  The  frequency  was  obtained  from  the  transducer  traces  of  the 
Individual  tubes  such  as  shown  in  Figure  21. 

The  value  of  the  quality  calculated  from  the  boiler  energy  balance  is 
considered  to  be  about  20  percent  lower  than  actually  existed.  The  reason 
for  this  disagreement  is  not  known,  since  the  heat  loss  calibration  was 
reproduced  five  times.  However,  the  fact  that  the  temperature  at  the  entrance 
to  the  cooler,  %diich  was  located  3  ft  and  two  90~degree  turns  beyond  the 
boiler  exit  (a  condition  which  tends  to  produce  equi 1 ibri urn),  remained  greater 
than  the  saturation  temperature  indicates  that  the  exit  quality  was  close  to 
100  percent. 

The  overall  coefficient  of  heat  transfer  was  calculated  from  the  net 
heat  input  to  the  Freon,  the  Freon  saturation  temperature,  and  the  measured 
air  temperatures.  The  heat  transfer  resistance  was  approximately  equally 
distributed  between  the  two  fluids. 

The  data  could  be  reproduced  without  difficulty,  indicating  that  all 
the  important  variables  were  under  control.  The  resulting  values  of  the 
controlling  dimensionless  parameters  were  calculated  from  the  measurements 
and  are  presented  in  Tables  5  and  6. 

The  data  for  each  run  corresponding  to  Figures  15  to  21  will  be  retained 
at  Ai Research  by  the  author. 
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FREON  BOILER  STABILITY  DATA  (FULL-LENGTH  TWISTED  TAPE) 


FREON  BOlLin  HOW  STABIlT'^v  DATA  {3/4-LENGTH  TWISTED  TAPE) 
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acquisttion  runs  were  completed,  10,000  ft  of  high-speed 
wtior.  pictures  were  obtained  of  both  stable  and  unstable  flow,  in  pfain^ 
tubes  and  tubes  with  twisted  tape  inserts.  The  motion  pictures  include 
equences  covering  a  wide  rarge  in  quality  and  showing  one  and  five  tubes. 

approximately  those  1000  ft  which  most 

w  it  be  and  to  produce  a  finished  motion  picture  which 

111  oe  available  on  load  from  the  author  in  approximately  one  year. 
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SECTION  8 


DISCUSSION  AND  CONCLUSIONS 


FLOW  STABiLirv  HEASURENSNTS 

The  measured  values  of  AP-j-p  range  from  0.2  to  1.6,  with  85  percent 

of  the  results  falling  below  a  value  of  1.0.  This  is  alnx}St  exactly  the  same 
range  of  this  parameter  that  was  obtained  when  Lowdermilk's  ( 10)  results  with 
water  were  presented  in  terms  of  this  parameter.  Although  the  range  oT  values 
of  the  other  dimensionless  parameters  in  Lowdermilk's  experiment  are  not 
available,  this  agreement  gives  added  confidence  for  the  use  of  Freon  flow 
stability  data  in  the  design  of  the  SNAP  SO/SPUR  boiler,  since  the  heat  of 
vaporization  of  water, like  that  of  potassium,  is  an  order  of  magnitude  greater 
than  that  of  Freon. 

To  clarify  the  relationship  between  the  liquid  pressure  drop  ratio  and 
t^n  other  dimensionless  parameters,  the  results  tabulated  in  Tables  5  and  6 
are  presented  graphically  in  Figures  24  to  27.  While  lines  have  been  drawn 
through  the  data  points  presented  In  the  figures,  the  location  and.  In  some 
cases,  even  the  direc«:lon  of  the  slope.  Is  very  uncertain,  since  only  three 
points  are  available  to  define  each  line.  Nonetheless,  certain  trends  seem 
to  exist. 

« 

Figures  24  tnd  25  reveal  that  the  liquid  orifice  pressure  drop  required 
to  produce  stable  flow  decreases  as  the  fractional  pressure  drop  Increases. 

This  Is  consistent  with  the  results  disclosed  by  the  literature  review,  and 
discussed  In  the  paragraphs  on  Inlet  Velocity  earlier  in  this  report,  which 
omclude  that  Increasing  Inlet  velocity  or  flow  rate  tends  to  stabilize  the 
system.  This  conclusion  Is  somewhat  weakened  by  the  reverse  slope  shown  in 
Flr**'e  25b  for  a  density  ratio  of  150,  and  by  the  scatter  In  Figure  25a  for 
the  same  density  ratio.  The  experimental  results  of  this  study,  however, 
when  combined  with  the  opinions  of  others,  leads  to  the  conclusion  that  an 
increase  in  fractional  pressure  drop  tends  to  stabilize  the  flow  in  a  forced- 
convection  boiler. 

Figures  26  and  27  show  a  tendency  for  ii  creased  stability  with  an 
Increase  in  subcooiing;  this  is  inconsistent  with  the  results  of  most  previous 
investigators.  Hoover,  Gouse  (14)  observed  in  his  recent  experiments  that 
there  was  a  range  of  subcooli'g  which  produced  instabili  y  and  that  increasing 
subcooling  beyond  this  range  tends  to  stabilize  the  system.  The  results  of 
the  present  experiment  tend  to  conflm  the  conclusion  that  increased  subcooling 
may  stabilize  as  well  as  destabilize. 

The  effect  of  density  ratio  (pressure  level;  is  best  displayed  in 
Figures  24  and  25.  For  the  range  of  density  ratio  investigated,  the  scatter 
in  the  data  almost  disguises  the  effect;  however,  increased  density  ratio 
appears  to  require  a  greater  orifice  pressure  drop.  This  conclusion  is 
considered  to  be  well  established  by  previous  work;  the  fact  that  the  present 
experimental  results  do  not  show  a  reverse  trend  is  considered  sufficient 
confirmation. 
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24.  Effect  of  Fractional  Pressure  Drop  with  3/4  Tape 


Effect  of  Fractional  Pressure  Drop  with  Full  Tape 


Figure  26.  Effect  of  Subcoot Ing  with  3/4  Length  Tape 
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igure  27.  Effect  of  Subcooting  with  Full  Length  Tape 


f  The  effect  of  exit  quality  was  not  intentionally  investiqated  in  this 

j  experiment.  However,  the  30-percent  variation  in  exit  quality  which  e^nsts 
I  in  the  results  is  probably  one  reason  for  the  scatter  in  the  data  exhibited 

I  in  the  figures,  since  the  accuracy  of  the  data  is  well  within  the  scatter 

i  that  is  evident. 


I  No  consistent  difference  is  discernible  between  the  results  for  the 

I  full-length  or  3/4-iength  twisted  tape  insert-  Any  effect  may  be  overcome 

I  by  variations  produced  by  different  exit  qualities  at  corresponding  points 

I  The  temperature  into  the  cooler,  shown  as  component  3  in  Figure  6,  is 

I  included  in  the  data  recorded  in  Tables  5  and  6  since  it  illustrates  the 

1  complexity  and  uncertainty  of  the  forced-convection  vaporization  process. 

I  This' temperature  generally  differed  from  the  boiler  exit  temperatuie,  which 
I  was  typically  30®F  greater  than  the  saturation  temperature;  in  some  ca‘^es, 
i  it  was  greater  and  in  other  cases  less.  All  of  these  observations  can  be 

I  explained  by  the  fact  that  the  vapor  in  high  quality  mist  flow  can  be  greatly 

I  superheated,  even  while  liquid  exis^  i  the  flow.  The  heat  first  is  trar^s- 

I  ferred  from  the  wall  to  the  vapor,*  len  vaporization  occurs  in  the  interior 

I  of  the  tube  due  to  heat  transfer  fr.  ^  ic  vapor  to  the  liquid.  Depending 

I  upon  the  location  of  the  thermocouple  and  the  way  in  which  the  superheated 

I  vapor  and  saturated  liquid  impinge  upon  the  thermocouple,  a  range  of  temper- 
I  atures  between  the  extreme  values  existing  in  the  tube  can  be  recorded.  Also, 

I  as  vaporization  occurs  in  the  unheated  length  between  the  boiler  and  the  cooler 

I  due  to  heat  transfer  between  the  vapor  and  the  liquid,  the  temperature  of  the 

I  vapor  tends  to  decrease.  The  multiplicity  of  nonequi 1 ibrium  temperatures 

t  which  it  is  possible  to  measure  in  high  quality,  two-phase  flow  gives  rise 
>  to  the  curious  combinations  of  data  recorded  in  these  experiments. 

I  PHOTOGRAPHIC  OBSERVATIONS 

I  The  high-speed  motion  pictures  provided  the  following  qualitative  inror. 

mation  about  two-phase  flow  in  a  tube  with  twisted  tape  inserts.  Beyond  the 
I  point  at  which  the  tube  wall  becomes  dry,  a  significant  quantity  of  liquid 
continues  to  flow  along  the  center  of  the  tape.  This  quantity  of  liquid 
j  appears  to  resist  the  tendency  to  centrifuge  to  the  tube  wall,  and  much  of 
I  it  persists  to  the  end  of  the  tape.  Interruptions  in  the  tape,  which  would 
q  Serve  to  shatter  this  liquid  film,  would  probably  cause  much  of  the  liquid 
in  the  film  to  centrifuge  to  the  wall  where  it  would  vaporize. 


The  photographs  also  confirmed  the  earlier  results  (with  plain  tubes) 
that,  even  when  the  flow  into  each  tube  was  steady,  the  exit  flow,  and  there¬ 
fore  quality,  always  fluctuated.  This  was  true  even  though  the  transducer 
traces  of  the  overall  boiler  pressure  drop  and  individual  tube  pressure  drop 
exhibited  no  fluctuations.  It  was  evident  from  the  visual  observations  that 
these  stable  flow,  exit  quality  oscillations  were  the  result  of  the  plug- 
:  annular  flow  regime  transition.  The  plugs,  which  exist  at  low  quality,  grow 
very  rapidly  as  they  flow  along  the  tube  and  additional  vapor  is  generated. 

Even  after  the  vapor  plugs  join, due  to  their  growth,  to  produce  annular  flow, 
a  disproportionate  amount  of  liquid  continues  to  exist  where  two  plugs  have 
P  joined.  This  liquid  is  pushed  along  the  tube  by  the  vapor,  and  probably 
r  much  of  it  ends  up  as  dispersed  droplets  in  the  vapor  core.  The  disproportionate 


56 


4. 


I 

I 

I 

I 

I 

1 

i 


I 

{ 


I 

I 

I 


amount  of  iiquid  which  exists  between  two  vapor  plugs  is  not  distributed 
uniformly  along  the  tube  after  annular  flow  occurs.  Consequently,  exit  flow 
and  quality  fluctuations  from  individual  tubes,  which  manifest  themselves 
visually  as  periodic  surges  of  mist  at  the  end  of  the  heated  length,  occur 
even  in  stable  two-phase  flow. 

The  distance  from  the  end  of  the  3/4- length  twisted  tape  to  the  end  of 
the  heated  length  was  approximately  10  in,  which  is  an  L/D  ratio  of  40.  In 
spite  of  this,  the  swirl  flow  produced  by  the  tape  persisted  to  the  end  of 
the  tube  in  all  runs.  There  .vas  a  visible  decrease  in  the  angle  between  the 
flow  direction  and  the  tube  axis  as  the  end  of  the  tube  was  reached,  but 
an  effective  amount  of  swirl  remained  at  the  end  of  the  tube. 

The  motion  pictures  provided  some  very  interesting  sequences  of  sub¬ 
cooled,  forced-convection  boiling.  Bubblers  nucleated  on  the  wall,  grew  to 
about  0.020  in.  in  diameter,  and  then  broke  away  from  the  wall  due  to  the 
inertia  caused  by  their  growth  rate.  After  breaking  away  from  the  wall  and 
penetrating  the  subcooled  liquid  core  in  the  interior  of  the  pipe,  the  bubbles 
very  rapidly  decreased  in  size  as  they  condensed,  until  only  a  smalt  spot 
was  visible  where  the  bubble  had  been  located.  When  this  small  spot  mcved 
downstream  and  reached  the  point  where  the  liquid  exceeded  the  saturation 
temperature,  it  again  grew  into  a  bubble  in  the  interior  of  the  tube.  Other 
forms  of  subcooled  boiling  were  observed  in  which  the  bubble  collapsed  while 
it  was  still  attached  to  the  wall,  but  the  process  described  above  was 
particularly  dramatic  and  interesting. 

SUHMARY  OF  CONCLUSIONS 

The  following  conclusions  appear  to  .Se  valid  for  the  conditions  investi¬ 
gated  in  this  experiment. 

1.  The  ratio  of  liquid  orifice  pressure  drop  to  two-phase  pressure 
drop  required  to  produce  stable  flow  for  SNAP"50/SPUR  boiler 
conditions  is  0.2  to  1.6.  To  be  conservative,  a  value  between 

1  and  2  should  be  used. 

2.  The  flow  tends  to  become  more  stable  as:  (a)  fractional 
pressure  drop  is  increased,  (b)  entrance  subcooling  is 
increased,  and  (c)  density  ratio  is  decreased. 

3.  The  vapor  in  mist  flow  Is  always  superheated,  which  makes  temperaf  re 
measurements  very  hard  to  interpret. 

4.  Exit  flow  and  quality  oscillations  from  individual  tubes  exist 
due  to  the  plug-annular  flow  regime  transition,  even  when  the 
flow  into  the  tube  is  steady. 

5.  Swirl  flow  persists  for  a  length  equal  to  more  than  40  diameters 
after  the  end  of  the  twisted  tape  in  high  qua1’'-y  two-phase  flow. 

6.  A  significant  quantity  of  iiquid  continues  to  flow  along  the 
center  of  a  plain,  twisted  tape  after  the  tube  wall  becomes  dry. 


57 


! 

I 


REFERENCES 


I 


1.  J.  G.  Collier,  "A  Review  of  Two-Phase  Heat  Transfer  1935-1957,"  AERE 
Report  CE/R  2496,  Harwell,  Berkshire,  England,  1958. 

2.  Peter  Griffith,  "Two-Phase  Flow  in  Pipes,"  Massachusetts  Institute  of 
Technology,  Course  2.529  Lecture  Notes,  July  1958. 

3.  Earl  t^uandt,  "Analysis  of  Gas-Liquid  Flow  Patterns,"  A.I.Ch.E.  Preprint 
No.  47,  Sixth  National  Heat  Transfer  Conference,  August  1963. 

4.  C.  E.  Oengler,  and  J.  N.  Addoms,  "Heat  Transfer  Mechanism  for  Vaporization 
of  Water  in  a  Vertical  Tube,"  Chemical  Engineering  Progress  Symposium 
Series,  No.  18,  Vol  52,  1956. 

5.  H.  Groothius  and  W.  P.  Hendal,  "Heat  Transfer  in  Tvx)-Phase  Flow," 

Chemical  Engineering  Science,  Vol  II,  pp  212  to  220,  1959. 

6.  P.  Sachs  and  R.A.K.  Long,  "A  Correlation  for  Heat  Transfer  in  Stratified 
Two-Phase  Flow  with  Vaporization,"  Int.  J.  of  Heat  and  Mass  Transfer, 

Vol  2,  pp  222-230,  1961. 

7.  John  C.  Chen,  "A  Correlation  for  Boiling  Heat  Transfer  to  Saturated 
Fluids  in  Convective  Flow,"  A.S.M.E.  Paper  63-HT-34,  August  1963. 

8.  P.  J.  Berenson  and  R.  A.  Stone,  "A  Photographic  Study  of  the  Meehans im 
of  Forced  Convection  Vaporization,"  A.I.Ch.E.,  Paper  No.  21,  Fifty-First 
National  Meeting,  1963,  to  be  published  in  I96'»  A.I.Ch.E.  Heat  Transfer 
Symposium. 

9.  J.  W.  Miles,  "The  Hydrodynamic  Stability  of  a  Thin  Film  of  Liquid  in 
Uniform  Shearing  Motion,"  J.  of  Fluid  Mechanics,  Vol  8,  1961,  pp  592-610- 

10.  W.  H.  Lowdermilk,  C-  0-  Lanzo  and  B-  L.  Siegel,  "Investigation  of  Boiling 
Burnout  and  Flow  Stability  for  Water  Flowing  In  Tubes,”  NACA  TN  4382, 

September  1958- 

11.  E.  S.  Beckjord,  "The  Stability  of  Two-Phase  Loops  and  Response  to 
Ships  Motion,"  GEAP  3493  (Rev.  l),  September  i960. 

12.  A.  P.  Fraas,  "Flow  Stability  in  Heat  Transfer  Matrices  Under  Boiling 
Conditions,"  Addendum  No.  J,  ORNL  59-11-1,  Oak  Ridge,  Tennessee  March  22,  1961 

i5.  W.  R.  Gonbill  and  N.  D.  Greene,  "Boiling  Burnout  with  Water  in  Vortex 
Flmi<, "  Chemical  Engineering  Progress,  Vol  54,  No-  lO,  October  1958 
pp  68-76 

14.  S.  W.  Gouse,  Jr.,  and  C.  C.  Andrysiak,  "Flow  Oscillations  in  a  Closed 
Loop  with  Transparent,  Parallel,  Vertical,  Heated  Channels,"  Report 
No.  8973-2,  Dept,  of  Mech.  Eng.,  Cambridge,  Mass.  June  1963. 


58 


REFERENCES  (Continued) 


15.  E.  R.  Quandt,  "Analysis  and  Heasurement  of  F’ow  Osci  Hat  ions, ’’  Chemicai 
Engineering  Progress  Symposiun  Series  No.  32,  Vol  57,  i96l,  pp  I! I- 126- 

16-  J.  E.  Meyer  and  R.  P.  Rose,  "AppJication  of  a  Momentum  Integral  Model 
to  the  Study  of  Parallel  Channel  Boiling  Flow  Oscillations,"  ASME 
Paper  No.  62-HT-41,  August  1962. 

17.  A.  L.  Blubaugh  and  E.  R.  Quandt,  "Analysis  and  Measurement  of  Flow 
Oscillations,"  (Westinghouse*Betlis  Atomic  Pow*er  Lab)  WAPD-AD-TH  538, 
p  72,  November  i960. 

i8-  0.  J-  Handler,  A.  S.  Rathbun,  N.  E.  Van  Huff,  and  A.  Weiss,  "Natural 

Circulation  Tests  with  Water  at  800  to  2000  psia  Under  Nonboiling, 

Local  Boiling,  and  Bulk  Boiling  Conditions,"  Journal  of  Heat  Transfer 
83  C  p  261,  August  1961. 

19.  A.  H.  Stenning,  "Instabilities  in  the  Flow  of  a  Boiling  Liquid,"  ASME 
Paper  62'WA-i55,  November  1962. 

20.  H.  R.  Wright  and  H.  F.  Poppendiek,  "Unsteady  State  Floi-^  in  Boiling 
Heat  Transfer  Systems,"  Technical  Progress  Report  No.  2,  GLM-15, 
Geoscience  Ltd.,  Solana  Beach,  Calif.,  October  1962. 

21.  A.  B.  Jones,  "Hydrodynamic  Stability  of  a  Boiling  Channel,"  Knolls 
Atomic  Power  Laboratory  Report  KAPL  -  2178,  October  1961. 

22.  G.  B.  Wallis  and  J.  H.  Heasley,  "Oscillations  in  Two-Phase  Flovj," 

The  American  Society  of  Mechanical  Engineers,  83,  C.  No.  3,  August  1961. 

23.  W.  H.  Lowdermilk,  C.  0.  Lanzo  and  B.  L.  Siegel,  "Investigation  of 
Boiling  Burnout  and  Flow  Stability  For  Water  Flowing  in  Tubes,"  NACA 
TN  4382,  September  1958. 

24.  S.  Octrach,  and  A.  Koester,  "Film  Instabilit;  es  in  Two-Phase  FIov/s," 
A.I.Ch.E.  Preprint  No.  45,  Sixth  National  Heat  Transfer  Conference, 
August  1963. 

25.  N.  Adorni,  I.  Casagrande,  L.  Cravarolo,  A.  Hassid,  and  H.  Silvestri, 
"Experimental  Data  on  Two-Phase  Adiabatic  Flow:  Liquid  Film  Thickness, 
Phase  and  Velocity  Distribution,  Pressure  Drops  in  Vertical  Gas-Liquid 
Flow,"  CISE  R35  (Euroatom),  March  1961. 

26.  H.  S-  Isbin,  R.  H.  Moen,  and  0.  R.  Mosher,  "Two-Phase  Pressure  Drops," 
Univ.  of  Minn.,  Tech.  Info.  Services,  Oak  Ridge,  Tenn.,  AECU  2994, 
November  1954. 

27.  H.  S.  Isbin,  "Comments  on  a  Paper  by  Hendler,  et  al.,"  Journal  of  Heat 
Transfer,  83  C,  p  261,  August  1961- 


59 


REFERENCES  (Continued) 


28.  R.  Delagre  and  J.  P.  Saunier.  "O’Acoulement  Du  Fluide  Dans  Un  Canal 
Chauffant  RefroidI  a  L*Eau — R,”  Commissar i at  A  L'Energie  Atomique, 

Centre  0  Etudes  Nudeaires^  Saclay,  France,  pp  55  (i96l). 

29.  A.  E.  Bergies  and  W.  M.  Rohsenow,  "Forced  Convection  Surface  Boiling 
Heat  Transfer  and  Burnout  in  Tubes  of  Small  Diameter,"  Tech.  Report 
No.  8767-21,  Dept,  of  Mech.  Eng.,  HIT 

50.  B.  Hatzner,  "Basic  Experimental  Studies  on  Boiling  Fluid  Flow  and  Heat 
Transfer  at  Elevated  Pressures,"  Columbia  Univ.,  Chem.  Eng.  MPR  1161, 

November  30,  1961. 

31.  S.  Stephen  Papelt,  "An  Instability  Effect  in  Two-Phase  Heat  Transfer  for 
Subcooled  Water  Flowing  Under  Conditions  of  Zero  F  *ty,"  ARS,  Space 
Power  Systems  Conference,  Santa  Monica,  Calif.,  Se,.  ..ember  25-28,  1962. 

32.  N.  A.  Radovcich  and  R.  Moissis,  "The  Transistion  from  Two-Phase  Bubble 
Flow  to  Slug  Flow,"  MIT  Tech.  Report  No.  7-7673-22,  June,  1962. 

33.  A.  P,  Fraas,  "Flow  Stability  in  Heat  Transfer  Matrices  Under  Boil'ng 
Conditions,"  ORNL  59-1 l-J,  Oak  Ridge,  Tennessee,  November  I,  1959 

34.  J.  G.  Collier  and  D.  J.  Pulling,  "Instabilities  in  Two-Phase  Flow," 

A.E.R.E.  Report  Ml  105,  Atomic  Energy  Research  Establishment,  Harwell, 
Berkshire,  England,  1963. 

35.  H  F.  Poppendiek,  and  C.  L.  Rickard,  "Remarks  on  Transient  Pressures 
Within  the  Core  of  the  Research  Reactor  Under  Power  Excursions," 

GAMD  176,  General  Atomics,  San  Diego,  Calif.,  July  1957. 

36.  E.  H.  Wissler,  H.  S.  Isbin,  and  N.  R.  Amundson,  "Oscillatory  Behavior 
of  a  Two-Phase  Natural  Circulation  Loop,"  The  American  Institute  of 
Chemical  Engineer’s  Journal  No.  2,  pp  157-162,  1956. 

37.  A.  R.  Gruber,  and  S.  C.  Hyman,  "Flow  Distribution  Among  Parallel  Heated 
Channels,"  The  American  Institute  of  Chemical  Engineer’s  Journal  No.  2, 
p  199,  1956. 

38.  F.  Schroeder,  S.  G.  Forbes,  W.  E.  Nyer,  F.  L.  Bentzen,  and  G.  0.  Bright, 
"Experimental  Study  of  Transient  Behavior  in  a  Subcooled  Water-Moderated 
Reactor,"  Nuclear  Science  and  Engineering  2,  No.  !,  pp  96-115,  February  1957. 

39.  W.  H.  McAdams,  Heat  Transmission,  McGraw-Hill  Book  Co.,  Inc.,  New  York 
1954,  Chapter  5. 


60 


appendix 


s 


II 


A  PHOTOGiMPHlC  STUDY  OT  THE  MECHANISM 
OF  FOACEO  CONVECTION  VAPORIZATION 


by 

P*  J*  Berenson 
R«  A*  Stone 

AtResearch  Manufacturing  Company 
Los  Angeles  9,  California 


ABSTRACT 

The  vaporizatic-i  of  Freon  113  flowing  Inside  a  horizontal  pyrex  tube 
was  studied  with  the  aid  of  7000  frame  per  second  motion  pictures.  Heat 
fluxes  up  to  40^000  Btu/hr-ft*  were  obtained  using  800*F  air  flowing  in  the 
space  between  the  glass  tube  and  a  rectangular  quartz  duct  as  the  heat  source. 
The  Freon  entered  the  tube  subcooled  and  the  exit  condition  varied  from 
inderate  quel  I ty  to  superheated  vapor.  The  two  (rfiase  flow  patterns  and  heat 
transfer  leechantsfes  revealed  In  the  fdiotographs  are  discussed.  The  photographs 
lead  to  the  conclusion  that  the  annular-mist  flow  transition  (forced  convection 
vaporization  burnout  condition)  results  simply  from  the  vaporization  of  the 
liquid  annulus  while  liquid  droplets  remain  dispersed  throughout  the  vapor 
core. 


1  INTRODUCTION 

The  ability  to  accurately  predict  the  heat  transfer  and  pressure  drop 
characteristics  of  the  forced  convection  vaporization  process  has  become  in¬ 
creasingly  important  in  recent  years.  Boiling  nuclear  reactors^  Rankine 
cycle  space  power  systems^  environmental  control  systems  for  spacecraft#  and 
other  advanced  systems  depend  upon  this  process. 

A  large  number  of  both  experimental  and  analytical  investigations  of 
forced  convection  vaporization  have  been  perforim.’*  during  the  past  twenty 
years.  Reference  I  contains  an  excellent  review  of  the  state-of-the-art  up 
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to  1957.  Considering  the  amount  of  ¥K>rk  that  has  been  done>  relatively  little 
progress  has  been  made  toward  predicting  the  characteristics  of  forced  convec* 
tion  vaporization.  This  is  due  to  the  great  number  of  variables  upon  which 
the  process  depends,  compounded  by  the  complexity  of  the  various  t%tfO*phase 
flow  patterns  which  occur  as  the  quality  increases  along  the  tube. 

Because  of  the  large  differences  in  the  various  flow  patterns  which  occur 
as  the  volume  fraction  of  the  vapor  increases  along  the  tube,  it  is  unlikely 
that  a  single  correlation  can  adequately  predict  the  heat  transfer  characteris¬ 
tics  of  forced  convection  vaporization  from  0  to  100  percent  quality.  It  is 
probable  that  a  separate  correlation  will  have  to  be  developed  for  each 
flow  pattern.  To  do  this  requires  a  familiarity  with  the  characteristics 
of  each  flow  pattern  so  that  an  analytical  model  can  be  formulated  to  yield 
correlations  for  the  heat  transfer  coefficient  in  each  regime  and  the  quality 
at  which  flow  regime  transitions  occur. 

Visual  studies  of  adiabatic,  two-phase,  tt«o-cofflponent  flow  have  provided 
information  on  the  flow  patterns.  Attempts  have  been  made  to  develop  flow 
regime  maps  which  t«ouid  make  it  possible  to  predict  the  transition  conditions 
(2).  While  the  ^iow  regimes  depend  to  some  extent  on  whether  the  flow  is 
vertical  or  horizontal,  for  both  orientations  the  flow  regimes  shown  in  Figure? 
appear  to  be  that  most  important  for  relatively  high  flow  rates.  The  flow 
patterns  shown  are  as  follows: 

1.  Bubbly  Flow — This  pattern  occurs  at  very  low  quality  and  consists 
of  individual  bubbles  of  vapor  entrained  in  the  liquid  flow.  The 
bubbles  are  formed  at  nucleation  sites  on  the  wall. 

2.  Plug  or  Slug  Flow — As  the  quality  Increases,  the  individual  bubbles 
agglomerate  to  form  plugs  or  slugs  of  vapor  which  periodically  pass 
a  given  point  on  the 

3.  Annular  Flow — At  higher  quality,  the  vapor  flows  in  a  high  velocity 
core  in  the  center  of  the  pipe  and  a  thin  film  of  liquid  covers 
the  pipe  wall.  Liquid  droplets  are  generally  dispersed  in  the 
vapor  core. 

Mist  Flow — At  high  quality,  the  liquid  film  on  the  wall  disappears 
and  all  the  liquid  is  dispersed  throughout  the  vapor  as  droplets. 
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Forced  convection  vapor jzst son  h«st  transfer  data  have  shuim  that  tne 
heat  transfer  coefficient  increases  as  the  quality  increases  (3«  4).  At  a 
critical  quality^  which  varies  considerably  with  fluid,  geonnetry,  and  heat 
flux,  there  is  a  sharp  decrease  in  the  heat  transfer  coefficient  toward  the 
value  characteristic  of  vapor  flow,  as  shown  in  Figure29.  When  the  quality 
reaches  100  percent,  the  coefficient  equals  the  vapor  heat  transfer  coefficient< 


It  is  often  of  great  importance  to  be  able  to  predict  the  quality  at 
which  the  sharp  decrease  in  heat  transfer  coefficient  occurs,  for,  if  the 
heat  flux  is  the  imposed  condition,  burnout  occurs  beyond  this  point*  It 
has  been  suggested  (3)  that  this  quality  corresponds  to  the  point  at  which 
the  transition  from  annular  flow  to  mist  flow  occurs.  This  is  very  likely, 
since  this  is  a  transition  from  a  I iquid*covered  wall  io  a  gas "blanketed 
wall,  which  is  similar  to  the  pool  boiling  burnout  condition  that  occurs  in 
changing  from  nucleate  to  film  boiling.  There  have  been  many  suggestions  as 
to  the  reason  for  the  transition  from  annular  to  mist  flow,  and  numerous 
correlations  of  the  forced  convection  burnout  condition  have  been  attempted. 

No  general  correlation  is  presently  available,  and  there  is  considerable 
confusion  regarding  the  transition  mechanism. 

Most  of  the  visual  studies  of  forced  convection  vaporization  have  been 
performed  under  adiabatic  conditions  because  of  the  relative  ease  of  this 
type  of  experiment.  In  addition,  most  of  the  heat  transfer  tests  have  used 
electrical  heating,  thereby  imposing  the  heat  flux  as  the  boundary  condition 
(5).  The  purpose  of  the  experiment  described  below  was  to  study  the  forced 
convection  vaporization  process  with  high-speed  motion  pictures  when  the 
overall  temperature  difference  is  the  imposed  boundary  condition.  In  addition. 
It  was  desired  to  obtain  motion  pictures  of  all  the  significant  flow  s>atierns, 
since  insufficient  emphasis  has  previously  been  placed  on  annular  and 
flow. 


The  motion  picture  supplement  to  this  paper,  showing  the  various  flow 
regimes,  is  available  on  loan  from  the  American  Institute  of  Chemical  Engineers 
and  t^  AiResearch  Manufacturing  Company. 
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11  APPARATUS  AND  PROCEDURE 


The  schematic  of  the  test  setup  isshown  in  Pigtire30and  the  test  section 
geometry  in  Figure  il-  The  test  section  consists  of  an  inner  Pyrex  tube  6  own 
10/  8  mm  00/  and  40  inches  long/  in  t^ir.h  Frecn  flows.  Air  flows  outside 
this  tube/  inside  a  20>inch-iong  rectangular  quartz  duct  whose  inside  ditnen- 
siops  are  3/8  inch  by  1-3/4  inch.  The  air  flow  Is  counter  to  the  Freon  flow. 

The  heated  length  of  the  test  section  is  approximate! y  24  inches/  consisting 
of  the  20  inches  of  quartz  tube  plus  approximately  2  inches  on  each  end 
inside  the  flange  connection.  The  Freon  tube  is  insulated  from  contact  with 
the  air  until  it  reaches  the  24-inch  test  section.  The  air  is  heated  in  a 
gas  furnace  to  approximately  800^F  and  metered  in  an  insulated  section  before 
it  passes  into  the  test  section. 

Manometer  pressure  taps  were  located  at  the  flanges  of  a  1-1/2  inch 
orifice  inside  a  2-inch  pipe  on  the  inlet  air  line  and  on  both  sides  of  the 
test  section.  These  permitted  measurement  of  the  air  flow  rate  through  the 
orifice  and  later  provided  a  means  of  estimating  the  air  leakage  from  the 
test  section.  Air  inlet  and  outlet  temperatures  were  measured  by  means  of 
thermocouples  in  the  air  stream,  freon  inlet  and  outlet  temperatures  were 
measured  by  thermocouples  located  in  mixing  sections  inside  the  Freon  tube. 

The  Freon  inlet  line  was  instrumented  with  an  expansion  and  contraction 
chamber/  with  transducer  taps  on  the  chamber  inlet  and  cutlet  lines  and  also 
in  the  central  stagnant  chamber.  Forward  flow  rate  was  indicated  as  an 
entrance  loss  in  the  latter  two  taps  in  the  expansion  chamber,  while  reverse 
flow  was  indicated  by  entrance  pressure  loss  in  the  first  two  taps  of  the 
test  section.  The  forward  and  reverse  pressure  drops  across  the  Freon  inlet 
flow  section  were  recorded  on  a  Sanborn  recorder.  Freon  flow  rate  was 
measured  by  using  a  stop  watch  to  measure  the  time  required  for  collection 
of  one  pound  of  Freon  from  the  condenser. 

A  WF-3  Fastax  camera  with  a  6-inch  f4.5  lens  was  mounted  approximately 
8  inches  from  the  quartz  tube,  shooting  horizontally  at  the  test  section  through 
a  6-inch  Aero  Ektar  supplementary  lens.  Ektachrome  ERB  (ASA  50)  film  was  used 
with  the  shutter  of  the  camera  wide  open.  At  this  lens  opening,  the  depth  of 
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field  over  which  the  cenera  was  In  focus  wa*i  approximately  one  mllliintter* 

The  test  section  was  backlighted  to  40,000  ft-candles  by  two  lOOO-watt 
guartz  Iodine  lamps  reflecting  from  white  cardboard*  The  film  speed  was 
approximately  7000  pictures  per  second;  each  motion  picture  run  took 
approximately  l-*/2  seconds,  and  approximately  half  of  this  time  was  required 
for  the  earner;  to  attain  full  speed. 

The  first  motion  pictures  obtained  showed  poor  resolution  between  Freon 
liquid  and  vapor.  As  a  resultj  Trace,  a  bright  red  Freon '■col  or  lug  agent 
used  in  detecting  leaks  in  Freon  air  conditioning  systems,  was  added  to  the 
Freon  supply  in  a  concentration  of  approximately  one  weight  percent. 

Subsequent  pictures  showed  excellent  resolMtion  and  detail  of  the  vaporiza¬ 
tion  process. 

Meat  input  to  the  Freon  was  determined  by  measuring  the  heat  lost  by 
the  air  when  no  Freon  was  flowing,  and  subtracting  that  quantity  from  the 
total  enthalpy  change  of  the  air  during  a  run.  It  was  later  learned  that 
there  was  a  leak  of  hot  air  through  the  flange  connections,  varying  from 
about  15  percent  to  25  percent  of  the  total  air  flow.  The  steady-state. heat 
loss  from  the  test  section  was  approximately  equal  to  the  heat  input  Into 
the  Freon;  as  a  result,  the  quality  at  the  test  section  was  later  estimated 
to  be  accurate  only  to  ±30  percent. 

At  the  start  of  each  run,  the  air  fl<xrf  through  the  unit  was  adjusted 
by  means  of  a  valve  in  the  inlet  air  line  to  the  desired  value.  The  Freon 
tank  was  then  pressurized  to  40  psig,  and  the  position  of  the  Freon  throttle 
valve  was  adjusted  until  the  proper  flow  rate  was  observed  leaving  the  conden¬ 
ser.  The  following  quantities  were  recorded:  air  Inlet  temperature,  air 
outlet  temperature.  Freon  inlet  temperature.  Freon  outlet  temperature,  test 
section  air  static  pressure,  test  section  air  /iP,  air  line  orifice  static 
pressure^  air  line  orifice  ZiP,  and  Freon  flow  rate  (the  time  it  took  for  I 
pound  of  Freon  to  be  collected  at  the  condenser).  Motion  pictures  were 
taken  at  the  same  spot  throughout;  approximately  6  inches  from  the  Frwn 
Inlet  to  the  heated  section  and  approximately  14  inches  from  the  pyrex  tube 
ent ranee. 
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III  DISCUSSION  OF  RESULTS 


A.  FLOW  STABILITY 

For  low  values  pressure  drop  across  the  throttle  valve  at  the  entrance 
to  the  heated  tube^  flow  instability  {i<e«f  osciilatvons  of  the  orifice 
pressure  drop)  was  observed*  This  is  consistent  with  the  results  of  other 
investigators  (6)*  Figure  32is  a  typical  pressure  transducer  trace  under 
these  conditions.  The  magnitude  of  the  pressure  and^  therefore,  flow  fluctua¬ 
tions  decreased  as  the  pressure  drop  across  the  inlet  throttle  valve  was 
increased  as  shown  in  Figure32.  The  flow  rate  fluctuations  were  apparently 
due  to  the  interaction  of  the  two-phase-flow  pressure  drop  characteristics 
with  the  compressible  volumes  in  the  system  (6).  When  the  inlet  throttle 
valve  is  choked,  the  flow  downstream  of  the  throttle  valve  is  insensitive  to 
conditions  upstream;  the  presence  of  the  compressible  volume  provided  by  the 
pressurizing  gas  on  top  of  the  Freon  no  longer  has  any  effect.  It  is  of 
interest  to  note  that  the  elimination  of  some  flexible  connections  between 
the  throttle  valve  and  the  heated  tube,  which  were  observed  to  be  periodically 
expanding,  yielded  increased  flow  stability.  At  all  heat  fluxes  and  flow 
rates  investigated,  the  inlet  flow  rate  could  be  made  steady  with  a  great 
enough  pressure  drop  across  the  inlet  throttle  valve*  Even  though  there 
was  steady  flow  into  the  tube,  the  flow  rate  and  temperature  at  the  exit 
fluctuated.  This  is  apparently  due  to  the  periodic  nature  of  plug  flow,  the 
compressible  volume  of  the  plugs,  and  the  fact  that  thermal  equilibrium  does 
not  exist  In  mist  flow. 

B.  FLOW  REGIMES 

A  summary  of  the  characteristics  of  the  flow  regimes  shown  In  the 
motion  picture  and  discussed  below  is  presented  in  Table  7. 
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TA8LE  7 

MOTION  PICTURE  SUMMARY 


Film  Strip 

1 

2 

3 

4 

5 

Flow  Regime 

Bubbly 

Plug 

Annular 

(Slug) 

Annular 

Annular*Nist 
Transit  ion 

I'reon  Flow  Rate, 
Ib/min  ±Z^ 

3.4 

1.0 

1.0 

1.0 

0.50 

Air  Flow  Rate, 
Ib/mln 

3.5  1O.2 

3.5  ±0.2 

3.3  ±0.2 

8  ±2 

8  ±2 

Heat  Flux 

BTU/hr  ft*  ±\0i 

20,000 

20,000 

20,000 

27,000 

35,000 

duality 

i30)( 

60*  subcooled 

a 

30|( 

I.  Bubbly  Flow 

A  photograph  typical  of  the  bubbly  floM  regime  shown  in  the  first  film 
strip  is  presented  in  Figure  3d  In  order  for  this  flow  pattern  to  exist,  the 
vapor  volume  fraction  must  be  low.  Because  of  the  high  ratio  of  IIquid*to* 
vapor  density  for  most  fluids  of  interest,  the  corresponding  mass  fraction 
(quality)  in  the  bubbly  flow  regime  is  extremely  small,  for  example,  much 
less  than  I  percent.  The  bubbly  flow  regime  can  occur  in  a  subcooled  liquid 
since  the  bubbles  nucleating  at  the  wall  may  move  into  the  liquid  core  be* 
cause  of  the  inertia  caused  by  their  growth  rates.  Thermal  equilibrium  does 
hot  exist  as  there  is  insufficient  time  for  the  bubbles  to  condense  before 
the  surrounding  liquid  is  heated  to  the  saturation  temperature  as  It  flows 
along  the  tube. 

In  the  bubbly  flow  regime,  the  vapor  volume  fraction  is  so  small  that 
it  has  little  effect  on  the  liquid  velocity.  Bubbles  grow  and  depart  from 
nucleation  sites  on  the  wall  the  same  as  in  nucleate  pool  boiling.  Depending 
upon  the  viall  and  fluid  temperatures  and  the  nucleation  charactet  !c$  jf 
the  wall,  the  nucleate  pool  boiling  heat  transfer  mechanism  may  or  may  not 
dominate  the  forced  convection  heat  transfer  to  the  liquid  occurring  bettween 
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the  bubble  nucieation  sites.  In  the  very  narrow  quality  range  over  which 
bubbly  flow  exists,  it  seems  reasonable  to  assume  that  the  heat  transfer 
can  be  predicted  by  superimposing  liquid  forced  convection  and  nucleate 
pool  boiling,  provided  the  superheat  ten^erature  difference  is  not  great 
enough  to  produce  film  boiling  (7). 

2.  Plug  Flow 

As  tne  vapor  volume  fraction  or  the  L/0  ratio  increases,  the  individual 
bubbles  begin  to  agglomerate  and  form  plugs  or  slugs  of  vapor,  as  shown  in 
Figure  34and  the  second  film  strip.  These  large  volumes  of  vapor  occur 
periodically  along  the  tube.  Even  though  the  volume  quality  has  increased 
to  approximately  lb  percent  in  the  film  strip,  the  corresponding  mass  fraction, 
or  quality,  is  approximately  0.1  percent.  The  presence  of  the  larger  volume 
fraction  of  vapor  begins  to  cause  a  significant  increase  in  the  liquid 
velocity  in  the  plug  flow  regime. 

The  plugs  of  vapor  are  compressible  volumes  which  make  possible  the  flow 
oscillations  within  the  tube  that  are  :>bserved  in  the  film  strip  even  though 
the  flow  entering  the  tube  is  steady.  Actual  flow  reversal  is  shown  in  the 
movie.  It  is  possible  that  the  periodic  nature  of  plug  flow,  combined  with 
the  compressible  volume  It  provides,  is  the  source  of  many  of  the  exit  flow 
oscillations  »d>ich  occur  in  forced  convection  vaporization  inasmuch  as  all 
other  flow  regimes  are  steady. 

The  film  strip  shows  that,  for  the  conditions  prevailing  in  this  expeii- 
ment,  bubbles  nucleate  and  depart  from  the  wall  at  numerous  sites  inside  the 
tube  in  the  plug  flow  regime.  The  heat  transfer  mechanism  in  the  plug  flow 
regime  is  probably  the  same  as  in  the  bubbly  flow  regime:  a  superposition 
of  forced  convection  to  a  liquid  and  nucleate  pool  boiling. 

While  both  the  bubbly  and  plug  flow  regimes  are  very  interesting,  it 
must  be  emphasized  that,  for  the  density  ratios  of  genera!  interest  in  forced 
convection  evaporators,  they  occur  up  to  only  very  low  quality,  so  that,  if 
one  is  vaporizing  to  high  quality,  these  flow  regir.es  have  a  relatively  small 
effect  on  the  heat  transfer  process.  They  only  become  important  if  the  tem¬ 
perature  difference  is  great  enough  to  cause  film  boiling  or  because  of  the 
flow  oscillations  produced  in  the  plug  flow  regime. 
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5«  Autiuiaf  Flow 

The  two  film  strips  showing  annular  flow  were  taken  without  moving  the 
camera.  A  typical  photograph  of  annular  flow  is  presented  in  Figure  5j.  The 
waves  present  on  the  liquid-vapor  interface  just  inside  the  tube  can  be  de¬ 
tected,  but  the  presence  of  the  continuous  two*phase  interface  makes  it  im¬ 
possible  to  see  into  the  vapor  core.  The  first  film  strip  is  at  a  quality 
of  approximately  I  percent,  while  that  of  the  second  film  strip  is  approxi¬ 
mately  5  percent.  The  increase  in  quality  was  produced  by  maintaining  the 
Freon  flow  and  increasing  the  hot  air  flow.  The  dark  spots  on  the  tube  wall 
are  dried  dye  left  from  previous  tests;  they  have  no  effect  on  the  flow  or 
heat  transfer  process. 

In  the  first  film  strip,  some  of  the  liquid  stratifies  into  the  bottom 
of  the  tube.  Around  the  balance  of  the  tube,  a  thin  film  of  liquid  covers 
the  wall;  heat  is  transferred  by  conduction  through  this  liquid  film.  Al¬ 
though  there  are  a  number  of  active  bubble  nucleation  sites,  these  probably 
have  a  negligible  effect  on  the  heat  transfer  rate.  The  heat  transfer  in 
annular  flow  is  very  similar  to  that  for  condensation  inside  a  tube;  it  is 
probably  valid  to  neglect  the  heat  transfer  to  the  stratified  liquid  and 
assume  that  all  the  heat  is  transferred  by  conduction  through  the  thin 
1 iquid  annulus. 

Periodic  flow  surges  occur  in  this  film  strip  which  must  have  been 
initiated  in  the  plug  flow  regime,  since  the  flow  entering  the  tube  was 
steady.  The  flow  surges  are  detected  by  the  deepening  color,  which  indi¬ 
cates  an  increased  liquid  volume  fraction.  This  flow  regime  is  sometimes 
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called  slug  flow,  although  the  heat  transfer  mechanism  is  typical  of 
annular  flow. 

In  the  second  annular  flow  film  strip,  the  quality  is  approximately 


i 


5  percent  and  the  stratification  of  the  liquid  has  essentially  disappeared 
because  the  increased  inertia  forces  now  dominate  gravity.  The  camera 
location  was  not  changed  as  the  quality  was  increased  in  order  to  observe 
the  same  active  nucleation  site  which  existed  in  the  center  of  the  previous 
f i Im  strip. 
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Aithcugh  the  vapor  core  cannot  be  seen  through  the  liquid-vapor  interface 
which  forms  the  inner  boundary  of  the  liquid  annulus^  the  high  velocity  vapor 
core  contains  many  liquid  droplets.  The  liquid-vapor  interface  exists  just 
inside  the  tube,  where  the  waves  caused  by  the  vapor  shear  are  apparent. 

Once  ago  n,  while  some  nucleation  sites  exist,  the  dominant  heat  transfer 
mechanism  is  probably  conduction  through  the  thin  liquid  film.  The  vapor  is 
generated  primarily  by  vaporization  from  the  liquid-vapor  interface  inside 
the  tube  and  not  by  the  formation  of  bubbles  in  the  liquid.  There  is  a 
significant  amount  of  liquid  dispersed  throughout  the  vapor  core  as  droplets, 
in  addition  to  the  liquid  in  the  annulus  on  the  wail. 

The  flow  surges  stili  occur  at  this  quality.  While  the  charge  in  color 
of  the  fluid  during  the  surge  is  not  as- obvious  as  at  lower  quality,  the 
deepening  pink  indicates  that  the  liquid  volume  fraction  is  increasing. 

4.  Mist  Flow  (Forced  Convection  Vaporization  Burnout  Condition) 

The  observation  of  the  transition  from  annular  to  mist  flow  is  of  great 
interest,  since  this  is  presumably  the  point  at  which  the  evaporator  heat 
transfer  coefficient  experiences  the  sharp  decrease  observed.  Therefore, 
this  transition  is  the  cause  of  the  forced  convection  vaporization  burnout 
condition.  A  photograph  of  the  partially  dry  wall  taken  from  the  last  film 
strip  is  shown  in  Figure  3u. 

The  flow  surges  which  were  apparent  In  the  annular  flow  regime  continue 
to  persist  in  this  transition  regime.  That  is,  the  wall  in  front  of  the 
camera  alternately  becomes  wet  with  liquid  (along  the  bottom  half  in  the 
movie)  and  then  dry.  This  occurs  about  three  times  in  the  film  strip. 

The  red  line  about  halfway  up  the  tube  diameter  is  a  thin  strip  of  dye 
which  presumably  is  less  volatile  than  the  Freon  and  tends  to  collect  at 
high  quality  in  the  manner  shown.  Liquid  droplets  dispersed  throughout  the 
vapor  core  can  be  seen  where  the  wall  is  dry.  These  droplets  were  formed 
downstream  at  lower  quality  and  are  carried  along  by  the  vapor.  The  velocity 
of  these  droplets  is  much  higher  than  the  velocity  of  the  liquid  film,  and 
therefore  gravity  has  a  much  smaller  effect  on  them. 
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The  following  observations  of  the  wall  drying  process  are  significants  f 

i 

A  snMil  dry  spot  forms  and  grows  in  all  directions  as  the  liquid  vaporizes  f 

because  of  conduction  heat  transfer  through  the  liquid.  The  small  strips  | 

of  liquid  on  the  wall  remain  almost  stationary  relative  to  the  high  velocity  | 

vapor  and  the  liquid  droplets  in  the  vapor  core  which  can  be  seen  flowing  in  f 

I 

the  tube.  The  dominant  heat  transfer  mechanisn  is  conduction  through  the  I 

I 

liquid  film;  nucleation  may  produce  the  initial  dry  spot  on  the  wall,  but  ! 

has  a  small  effect  on  heat  transfer.  This  is  the  same  drying  process  that 
would  occur  with  a  thin  film  of  liquid  in  a  hot  pan  whose  temperature  is  not 
great  enough  to  cause  nucleate  boiling.  This  transition  occurs  three  times  f 

I 

in  the  film  stripy  as  described. 

Beyond  the  point  at  which  these  film  strips  were  toivcn,  t!^e  wall  is  > 

dry  except  when  the  liquid  droplets  in  the  vapor  core  impinge  upon  the  wall. 

Most  of  the  heat  transfer  is  from  the  wall  to  the  vapor.  After  the  heat  is 
transferred  to  the  vapor,  it  is  then  transferred  to  the  liquid  droplets;  the  I 

mist  flow  vaporization  process  actually  takes  place  in  the  interior  of  the 
pipe,  not  at  the  wall.  For  this  reason,  the  temperature  of  the  vapor  in 
the  mist  flow  regime  can  be  greater  than  the  saturation  temperature;  thermal 
equilibrium  may  not  exist  in  the  tube. 

The  above  description  of  the  transition  from  annular  to  mist  flow  is  in 
conflict  with  theories  which  postulate  that  the  transition  occurs  when  the 
thin  liquid  annulus  becomes  unstable  owing  to  the  effect  of  the  high  velocity 
vapor  core  (4).  This  was  not  observed  in  the  high  speed  photographs.  As  a 
result  of  the  visual  observations  of  the  transition  process,  it  is  concluded 
that  t‘ie  transition  is  simply  due  to  the  liquid  film  on  the  wall  vaporizing 
as  a  result  of  conduction  heat  transfer  through  the  liquid,  while  the  vapor 
core  contains  a  dispersion  of  liquid  droplets.  While  the  volume  fraction 
of  these  droplets  is  very  small,  they  account  for  a  substantial  mass  fraction 
because  of  the  high  ratio  of  the  I iquid-to>vapor  density. 

This  explan.  ^ion  is  consistent  with  a  theoretical  stability  analysis 
of  a  liquid  film  (8),  whose  predictions  have  been  partially  verified  by 
experiment  (9).  This  analysis  predicts  that  a  liquid  film  is  stable  for 
sufficiently  small  Reynold  numbers,  independent  of  the  vapor  velocity. 
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Since  the  liquid  film  Reynolds  number  decreases  as  the  quality  increases  in 
a  forced  convection  evaporator^  the  liquid  annulus  will  be  stable  at  suffi¬ 
ciently  high  quality,  independent  of  the  value  of  the  vapor  velocity.  There¬ 
fore,  the  high  quality  burnout  condition  could  not  be  due  to  an  instability 
of  the  liquid  film. 

While  annular  two-phase  flow  is  sometimes  idealized  so  that  all  the 
liquid  is  in  the  annular  film  and  all  the  vapor  in  the  high  velocity  core, 
observations  have  shown  that  droplets  of  liquid  generally  exist  in  the  vapor 
core.  These  droplets  may  or  may  not  join  the  liquid  annulus,  depending  on 
the  magnitude  of  the  turbulent  vapor  fluctuations  and  the  transit  time  of 
the  droplets  through  the  evaporator.  This  is  consistent  with  the  stability 
analysis  which  predicts  a  high  Reynolds  number  liquid  film,  such  as  exists 
at  low  quality,  may  be  unstable. 

In  summary,  the  liquid  droplets  dispersed  throughout  the  vapor  core 
are  produced  as  a  result  of  a  film  instability  at  low  quality.  The  tran¬ 
sition  from  annular  ^o  mist  flow  occurs  when  the  liquid  annulus  vaporizes 
as  a  result  of  conduction  through  the  thin  liquid  film,  while  liquid  drop¬ 
lets  remain  dispersed  throughout  the  vapor  core.  The  quality  at  the  tran¬ 
sition  is  determined  by  the  mass  of  liquid  contained  in  the  droplets  when 
the  transition  occurs. 


IV  SUMMARY  OF  CONCLUSIONS 

The  following  conclusions  appear  to  be  valid  based  upon  the  observations 
of  the  forced  convection  vaporization  process. 

1.  Throttling  at  the  entrance  of  a  forced  convection  evaporator  tube 
yields  steady  flow  into  the  tube. 

2.  Even  though  the  flow  into  the  tube  is  steady,  the  exit  flow  can 
fluctuate.  This  is  probably  due  primarily  to  the  periodic  nature 
of  plug  flow  combined  with  the  compressible  volume  within  the  tube 
provided  by  the  plug. 
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3.  NucUatlon  of  bubbUf  on  tho  ¥<a!!  !«as  observed  stheoever  tb*  wilS 
Mas  wety  in  all  flow  regtnas. 

4.  It  Is  reasonable  to  correlate  bubbly  and  plug-stu9  flow  by 
super i»|>osln9  liquid  forced  convection  and  nucleate  pool 
boiling  heat  fluxes. 

3*  The  dominant  heat  transfer  aiechanisai  in  annular  flow  is  conduction 
through  the  liquid  film  on  the  wall.  The  vapor  fonaatton  process 
occurs  primarily  at  the  annulus  liquid-vapor  interface  within  the 
tube#  not  at  the  wall  within  the  liquid  annulus. 

6.  The  transition  from  annular  to  mist  flow  (forced  convection  vaport- 
]:ation  burnout)  is  simply  the  result  of  the  annular  liquid  film 
Vaporizing  while  liquid  droplets#  fonaed  at  low  quality#  remain 
dispersed  throughout  the  vapor  core. 

7.  Departure  from  thermal  equilibrium  is  Im^rtant  In  forced  con¬ 
vection  vaporization;  bubbly  flow  can  exist  in  a  subcooled  liquid 
and  the  vapor  in  mist  flow  can  be  superheated. 
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Figure  28.  Two-Phase  Flow  Regimes 
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Figure  29.  Characteristic  Forced  Convection 
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Figure  30.  Test  Apperetus  Schemetic 
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Figure3l.  Test  Section  Geometry 
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Figure 32.  Typical  Preisure  Oscillations 


Figure  33.  Bubbly  Flow 
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